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Measurements of the Af, — > p£~vi and A;, — ¥ A c £~i>i decay rates can be used to determine 
the magnitudes of the CKM matrix elements V u b and V c i ,, provided that the relevant hadronic form 
factors are known. Here we present a precise calculation of these form factors using lattice QCD with 
2+1 flavors of dynamical domain-wall fermions. The b and c quarks are implemented with relativistic 
heavy-quark actions, allowing us to work directly at the physical heavy-quark masses. The lattice 
computation is performed for six different pion masses and two different lattice spacings, using 
gauge-field configurations generated by the RBC and UKQCD collaborations. The b —> u and b —> c 
currents are renormalized with a mostly nonperturbative method. We extrapolate the form factor 
results to the physical pion mass and the continuum limit, parametrizing the g 2 -dependence using 
^-expansions. The form factors are presented in such a way as to enable the correlated propagation 
of both statistical and systematic uncertainties into derived quantities such as differential decay 
rates and asymmetries. Using these form factors, we present predictions for the A & — > p£~vi and 
A;, —> A c l~vt differential and integrated decay rates. Combined with experimental data, our results 
enable determinations of \V u b\, |I4f>|, and IWt/V^I with theory uncertainties of 4.4%, 2.2%, and 
4.9%, respectively. 


I. INTRODUCTION 

To date, all direct determinations of the CKM matrix element magnitudes \V u b\ and \V c b\ were performed using 
measurements of B meson semileptonic or leptonic decays at e + e _ colliders. For both \V u b\ and \V c b\, there are 
tensions between the most precise extractions from exclusive and inclusive semileptonic B decays. The 2014 Review 
of Particle Physics lists [1] 


IK&lexci. = (3.28 ± 0.29) x 10~ 3 , IKfelexci. = (39.5 ± 0.8) x 10~ 3 , 

IKblincl. = (4.41 ± 0.15±g;if) x 10-3, |U cb | incL = (42.2 ± 0.7) x 10-3. ! 

The exclusive results in Eq. (1) are from the decays B -+ nil' and B —> D*£v (where £ = e, p) and use hadronic form 
factors from lattice QCD [2, 3]. The discrepancy between the exclusive and inclusive results is a long-standing puzzle 
in flavor physics [4-6], and right-handed currents beyond the Standard Model have been considered as a possible 
explanation [7-10]. New lattice QCD calculations of the B —> 7r form factors published recently yield somewhat 
higher values of |T4f,| eX ci. = (3.72 + 0.16) x 10~ 3 [11] and |I4f,| eX ci. = (3.61 + 0.32) x 10 -3 [12], but the latest analysis of 
B -+ D*tv using lattice QCD gives IV/blexci. = (39.04 + 0.75) [13] and slightly increases the exclusive-inclusive tension. 
Moreover, the current experimental results for the ratios of the B -+ D^tv and B -+ D^tv {£ = e,p) branching 
fractions differ from the Standard-Model expectation with a combined significance of 3.4cr [14]. 

On the experimental front, new results are expected from the future Belle II detector at the SuperKEKB e + e _ 
collider, and in the near future also from LHCb at the Large Hadron Collider. The LHCb Collaboration is currently 
analyzing the ratio of branching fractions of the baryonic b —t u and b —> c decays Af, — t pp~D^ and Af, -+ A c p~D^, 
with the aim of determining \V u b/V c b\ for the first time at a hadron collider. These decays were chosen over the 
more conventional B -+ 7 rpD and B —> DpD decays because, with the LHCb detector, final states containing protons 
are easier to identify than final states with pions [15]. Note that the production rate of Af, baryons at the LHC is 
remarkably high, equal to approximately 1/2 times the production rate of B° mesons [16]. 

The extraction of \V u b\ and \V c b\ (or their ratio) from the measured A& -+ and Af, -+ Ac/z - ^ branching 

fractions requires knowledge of the form factors describing the Af, —> p and Af, —>■ A c matrix elements of the relevant 
b —Y u and b —> c currents in the weak effective Hamiltonian. These form factors have been studied using sum rules 
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and quark models [17-28]. Nonperturbative QCD calculations of the Af, —> p and Ah —► A c form factors can be 
performed using lattice gauge theory. The first lattice QCD calculation of Ah —tp form factors, published in Ref. [29], 
employed static b quarks (i.e., leading-order heavy-quark effective theory) to simplify the analysis. The static limit 
reduces the number of independent Ah —> p form factors to two [30-32], but introduces systematic uncertainties of 
order Aqcd /mt and |p'|/mh in the Ah —> p 9^ differential decay rate (where p' is the momentum of the proton in 
the Ah rest frame). Here we present a new lattice calculation which improves upon Ref. [29] by replacing the static b 
quarks by relativistic b quarks, eliminating this systematic uncertainty. In addition to the six form factors describing 
the hadronic part of the decay Ah —► p p~9^ in fully relativistic QCD, we also compute the six analogous form factors 
for Ah —> A c p~ 9^ (note that early lattice studies of Ah —» A c form factors in the quenched approximation can be 
found in Refs. [33, 34]). Preliminary results from the present work were shown in Ref. [35]. 

In Sec. II we provide the definitions of the form factors employed here. The lattice actions and parameters, as well 
as the matching of the b —» u and b —> c currents from the lattice renormalization scheme to the continuum MS scheme 
are discussed in Sec. III. This calculation is based on the same lattice gauge-field ensembles as Ref. [29]; the ensembles 
include 2+1 flavor of dynamical domain-wall fermions and were generated by the RBC and UKQCD Collaborations 
[36] . Section IV explains our method for extracting the form factors from ratios of three-point and two-point correlation 
functions and removing excited-state contamination by extrapolating to infinite source-sink separation. Our fits of 
the quark-mass, lattice-spacing, and momentum-dependence of the form factors are discussed in Sec. V. The form 
factors in the physical limit are presented in terms of ^-expansion [37] parameters and their correlation matrices. Two 
different sets of parameters, referred to as the “nominal parameters” and the “higher-order parameters” are given. 
The nominal parameters are used to obtain the central values and statistical uncertainties of the form factors (and 
of derived quantities), while the higher-order parameters are used to calculate systematic uncertainties. In Sec. VI 
we then present predictions for the Af, -+ p£~9p and Af, — > A c £~9p differential and integrated decay rates using our 
form factors. Combined with experimental data, our results for the A b —> p fi and A;, —> A c /uz+ decay rates in the 
high-g 2 region will allow determinations of \V u b\ and \V c b\ with theory uncertainties of 4.4% and 2 . 2 %, respectively. 


II. DEFINITIONS OF THE FORM FACTORS 


Allowing for possible right-handed currents beyond the Standard Model, the effective weak Hamiltonian for b —► 
q t~ 9p transitions (where q = u,c) can be written as 


= G ^ V q L b [(! + £ g)<n^ b - (! - <+ ) 97^75 &] ^7m( 1 - 7s)^ 


( 2 ) 


(in the Standard Model, e q = 0 and V q r f } = V q i ,). To calculate the differential decay rate and other observables, we 
therefore need the hadronic matrix elements of the vector and axial vector currents, q'y^b and < 77 ^ 756 . In the following, 
we denote the final-state baryon by X (X = p , A c ). Lorentz and discrete symmetries imply that the matrix elements 
(AT|g7 M b |Af,) and (V| <77^75 fe| A&) can each be decomposed into three form factors. In this work we primarily use a 
helicity-based definition of the Af, —► X form factors, which was introduced in Ref. [38] and is given by 
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( 4 ) 


In these expressions, q = p — p' is the four-momentum transfer (whereas q is the wore quark field), and s± is defined 
as 


S± = {m Ab ± m x ) 2 - q-. 


( 5 ) 
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The form factors with subscripts 0, +, _L describe the contractions of the above matrix elements with virtual polar¬ 
ization vectors e* that are, respectively, time-like, longitudinal, and transverse to q 1 ''. Consequently, this choice of 
form factors leads to particularly simple expressions for observables such as the differential decay rate. Moreover, this 
choice simplifies the extraction of the form factors from correlation functions and clarifies the spin-parity quantum 
numbers of poles outside the physical kinematic region 0 < q 2 < {rn\ b — m x ) 2 ■ 

An alternate definition of the form factors that can be found in the literature (see, e.g., Ref. [27]) is the following: 


(X(p',s')\qYb\k b (p)) = u x {p',s') 
{X(p', s , )I?7 m 75 b\A b (p)) = u x (p',s') 




, W) M 


m A b 


q„ + 


“Ai 
pA(„ 2\ 


} A {q , )lf _gm ir ,^ + !im qr 

mx mx 


ua b (p,s), 


75 ua b {p,s), 


( 6 ) 

(7) 


where = §(7 M 7" — 7^7^) and, as before, q = p — p '. This choice decomposes the matrix elements into form 
factors of the first and second class according to Weinberg’s classification [39]. The second-class form factors and 
fY would vanish in the limit mb = m c (for At —> A c ) or mb = m u (for At —> p) [40]. In the following, we will refer 
to the form factors defined in Eqs. (6), (7) as “Weinberg form factors”. The helicity form factors are related to the 
Weinberg form factors as follows: 


U(q 2 ) 

= fi ( q 2 ) + —, q -r/ 2 V (? 2 )» 

TOAf, {rriA b + mx) 

(8) 

f±(q 2 ) 

-/r+)+ mA ‘ +mx /t+), 

mA b 

(9) 

fo(q 2 ) 

-fi(q 2 )+ , n/3 F (9 2 )» 

m Ab \jn Ab - m x ) 

(10) 

g+(q 2 ) 

- fi(q 2 ) , x f 2 (q 2 ), 

m.A b (m A b - mx) 

(11) 

g±(q 2 ) 

= V) _m A6 -TO. A 
mA b 

(12) 

go(q 2 ) 

- fi(q 2 ) , Jz(q 2 )- 

m Ab (rriA b + m x ) 

(13) 


These relations also demonstrate the following endpoint constraints for the helicity form factors: 


/o(0) = /+ (0), (14) 

So(0) = 5+(0)) (15) 

max) 9+{Q max)) (15) 

where q^ ax = ( mA b — m x ) 2 - At intermediate stages of our analysis of the lattice QCD data, it is beneficial to work 
with both definitions of the form factors. However, we perform the chiral/continuum/kinematic extrapolations only 
in the helicity basis. 


III. LATTICE ACTIONS AND CURRENTS 


This calculation is based on lattice gauge field ensembles generated by the RBC and UKQCD collaborations [36] 
with the Iwasaki gauge action [41, 42] and 2+1 flavors of dynamical domain-wall fermions [43-45]. We implement the 
light (u or d) valence quarks with the same domain-wall action that was used in generating the ensembles. Our analysis 
uses six different combinations of light-quark masses and lattice spacings as shown in Table I. These parameters are 
identical to those used in the earlier calculation of A b —► p£v form factors in Ref. [29] . However, instead of the static 
Eichten-Hill action [46] employed in Ref. [29], we now use anisotropic clover actions for the heavy (c and b) quarks 
[47-50]. These actions have the form 


sq = o, 4 J2q 


m Q +7oV 0 - ^V„ 2) 






i= 1 


i= 1 


i,j =1 


Q, (17) 


T— 7 ( 2 ) 

where Q is the lattice charm or bottom quark field, and V), ; are first- and second-order covariant lattice derivatives, 
and /*)„, is a lattice expression for the field-strength tensor (all of which are defined as in Ref. [51]). By suitably tuning 
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Set 

P 

Ng x At x As 

am 5 

(sea) 

am s 

(sea) 

a (fm) 

( v al) 

am u,d 

ml val) (MeV) 

N meas 

C14 

2.13 

24 3 x 64 x 16 

1.8 

0.04 

0.005 

0.1119(17) 

0.001 

245(4) 

2672 

C24 

2.13 

24 3 x 64 x 16 

1.8 

0.04 

0.005 

0.1119(17) 

0.002 

270(4) 

2676 

C54 

2.13 

24 3 x 64 x 16 

1.8 

0.04 

0.005 

0.1119(17) 

0.005 

336(5) 

2782 

F23 

2.25 

32 3 x 64 x 16 

1.8 

0.03 

0.004 

0.0849(12) 

0.002 

227(3) 

1907 

F43 

2.25 

32 s x 64 x 16 

1.8 

0.03 

0.004 

0.0849(12) 

0.004 

295(4) 

1917 

F63 

2.25 

32 s x 64 x 16 

1.8 

0.03 

0.006 

0.0848(17) 

0.006 

352(7) 

2782 


TABLE I. Parameters of the lattice gauge field ensembles [36] and light-quark propagators [29, 53]. The three groups of data 
sets {C14, C24, C54}, {F23,F43}, and {F63} correspond to three different ensembles of lattice gauge fields: one with a “coarse” 
lattice spacing a ta 0.11 fm, and two with “fine” lattice spacings a sa 0.085 fm (we use the lattice spacing values determined in 
Ref. [54]). Within each group, the valence-quark masses used for the propagators differ, resulting in different “valence 

pion masses” mi va1 ’ 1 ; the number of light-quark propagators used in each data set is denoted as A meas . 


Parameter 

coarse 

fine 

(6) 

am 

8.45 

3.99 

£(*>) 

3.1 

1.93 

c (b) 

U E, B 

5.8 

3.57 

(c) 

am Q 

0.1214 

-0.0045 

£(c) 

1.2362 

1.1281 

c (c) 

C E 

1.6650 

1.5311 

c (c) 

C B 

1.8409 

1.6232 


TABLE II. Parameters of the bottom and charm quark actions [51, 52]. 


the parameters z/, eg, cb as functions of am.Q, heavy-quark discretization errors proportional to powers of cluiq can 
be removed to all orders. The remaining discretization errors are of order a 2 |p| 2 , where |p| is the typical magnitude 
of the spatial momentum of the heavy quark inside the hadron. As the continuum limit a —i 0 is approached, the 
values v = 1 and Ce = cb = csw corresponding to the standard clover-improved Wilson action are recovered. For the 
bottom quark, we use the parameters that were tuned nonperturbatively by the RBC and UKQCD collaborations [51] 
using the condition that the action reproduces the correct spin-averaged B s meson mass and relativistic dispersion 
relation, as well as the correct B* — B s hyperfine splitting. For the charm quarks, we use the parameters from 
Ref. [52], where am,Q and v were tuned nonperturbatively to obtain the correct spin-averaged charmonium mass and 
relativistic dispersion relation, while ce and cb were set to mean-field improved tree-level predictions. Note that after 
the parameters were tuned in this way, the calculated charmonium hyperfine splittings were also in agreement with 
experiment [52]. The values of all lreavy-quark action parameters used here are given in Table II. 

We use a mostly nonperturbative method [55, 56] to match the b —I q (q = u,c) vector and axial vector currents 
from the lattice scheme to the continuum MS scheme. The renormalized currents in the MS scheme are written in 
terms of the lattice quark and gluon fields as 


Vo 


Ao 

Vi 

Ai 


7(99) y(bb) 

^V PVo 

7(99) y(bb) 

AV P A o 

Z^zWpVi 

z^zf\ Ai 


97o & + 2 a (cy 0 <77o7i^ jb + Cy 0 9 V jjoljb 
97075 b + 2 a (c ^ 0 q^ololj^ jb + c Ao 9^y 

qjib + 2 a (cy. qji'YjV jb + Cy. qv jji'Yjb + dy. 9 V ib + dy. qty t b 
97i75& + 2a (c A . 97i7s7i ^jb + c A . gV j'Yi'ys'ijb + d A . q'y^ib + d A . 9 ^ 75 ^ 


(18) 

(19) 

( 20 ) 
( 21 ) 


where Zy^ and Zy ^ are the matching factors of the flavor-conserving temporal vector currents 9709 and bjob, which 
are computed nonperturbatively using charge conservation. These nonperturbative factors provide the bulk of the 
renormalization, resulting in a much improved convergence of perturbation theory for the residual matching factors 
Pv^Ap- Above, i denotes the spatial components (i = 1,2,3), and the repeated index j is summed from 1 to 3. The 
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Parameter 

b —> c, coarse 

b — ¥ c, fine 

b —¥ u, coarse 

b —► u, fine 

p v ° 

0.9798(20) 

0.9848(15) 

1.02658(69) 

1.01661(52) 

Pa° 

1.0193(15) 

1.0112(29) 

1.02658(69) 

1.01661(52) 

Pvi 

1.0184(38) 

1.0162(41) 

0.99723(25) 

0.99398(12) 

pAi 

0.9866(33) 

0.9896(26) 

0.99723(25) 

0.99398(12) 

r R 

CyO 

0.0258(13) 

0.02873(99) 

0.0558(63) 

0.0547(64) 

r R 

C A 0 

0.03500(93) 

0.03285(87) 

0.0558(63) 

0.0547(64) 

CyO 

-0.0183(32) 

-0.0135(22) 

-0.0099(99) 

-0.0095(95) 

C A o 

-0.0205(17) 

-0.0155(27) 

-0.0099(99) 

-0.0095(95) 

c R 

L 'V3 

0.03192(70) 

0.0305(10) 

0.0485(27) 

0.0480(30) 

c R 

C A3 

0.0221(31) 

0.0237(20) 

0.0485(27) 

0.0480(30) 

C L . 

0.0088(22) 

0.0027(35) 

-0.0033(33) 

-0.0020(20) 

C L . 

-0.0002(31) 

-0.0020(20) 

-0.0033(33) 

-0.0020(20) 

d R 

u y3 

-0.0055(12) 

-0.0067(17) 

-0.00079(79) 

-0.0012(12) 

jR 

R A3 

0.0060(69) 

0.0044(40) 

0.00079(79) 

0 .0012(12) 

dyj 

-0.0176(44) 

-0.0043(69) 

0.0018(18) 

0.00047(48) 

<4 

0.0134(69) 

0.0106(40) 

-0.0018(18) 

-0.00047(48) 


TABLE III. Perturbative renormalization and improvement coefficients. The uncertainties given here include estimates of 
the missing higher-loop corrections as well as uncertainties from the numerical evaluation of the one-loop integrals. They are 
explained in more detail in the main text. 


lattice currents containing the lattice derivatives 

b = [U^(x)b{x + A) - U\{x - p)b{x - £)] , (22) 

9 ^7/x = -t [<?(a; + p)Ul(x) - 9^ - - A)] , (23) 

are needed to remove (H(a)-discretization errors from the currents. Time derivatives have been eliminated in Eqs. (18- 
21) using the equations of motion. We have computed Pv ,A as well as all of the 0(a)-improvement coefficients Cy A , 
d v ' A to one loop in mean-field improved lattice perturbation theory using the automated framework PhySyHCAl 
[57, 58]. The results are given in Table III. The central values are the average of plaquette and Landau-gauge mean- 
field improved results with perturbative expansion in = a -1 ) [51]. The uncertainties are the maximum of 

i) the difference of the respective mean-field improved results, ii) the numerical integration error, and iii) a power¬ 
counting estimate. For consistency with earlier stages of this project a different power-counting estimate is used for 
the b -A u and b —> c cases. For a perturbative quantity h with tree-level result h and full one-loop result h W we 
use (h^/h^ — l) 2 /)/ 1 ) as an error estimate for b u and (A 1 ) — h^)(a^g(p = a -1 )/ n) as an estimate for b —» c. 
For the p factors (but not for the a) -improvement coefficients) this estimate tends to be less conservative than 
estimates that are used in similar work [11, 12, 59]. The estimates of the combined uncertainty from the p factors 
and 0(a)-improvement coefficients, see Figs. 11, 13 and Table XII, agree well with similar work [11, 12, 59]. 

The nonperturbative matching factors Zy' 11 and Zy b ■* are given in Table IV. The light-quark and bottom-quark 
Zy U ' and Zy bS> were computed by the RBC and UKQCD collaborations [36, 60]. We determined the charm-quark 
Zy ^ using the method of Ref. [60] , by computing the following ratio of D s meson correlation functions without and 
with insertion of the current Jq = cyoc: 


^2(D s (xo + t,z) D\(x 0 ,x)) 

Z 

Y^(Ds(xo+t,z) J 0 (x 0 +t',y) Dl(x 0 ,x)) 
y,z 


Here, we used the following interpolating field with the quantum numbers of the D s meson, 


(24) 


D s = S75C, 


( 25 ) 
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FIG. 1. Numerical results for the ratio R ( cc ) (t, t') at t/a = 20 for the C54 data set (left) and at t/a = 24 for the F43 data set 

z v 

(right). The horizontal lines indicate the extracted values of Zy C \ 


Parameter 

coarse 

fine 

y(bb) 
ty (cc) 

z v 

ry(uu) 

z v 

10.037(34) 

1.35725(23) 

0.71651(46) 

5.270(13) 

1.18321(14) 

0.74475(12) 


TABLE IV. Nonperturbative renormalization factors of the flavor-conserving temporal vector currents. For Zy U \ we use the 
results in the chiral limit from Ref. [36]. For Zy b \ we use the results obtained in Ref. [60] on the coarse = 0.005 and 

fine am^ e ^ = 0.004 ensembles. 

where the tilde indicates gauge-covariant Gaussian smearing to suppress excited-state contamination. For large 
Euclidean time separations t, t' , and 1 1 — t'\, the ratio (24) becomes equal to Zy C \ Our numerical results for 
R„(cc)(t,t'), along with fits in the plateau region giving z[/ c \ are shown in Fig. 1 for the ensembles used in this 
calculation. 


IV. EXTRACTION OF THE FORM FACTORS FROM CORRELATION FUNCTIONS 

In this section we explain how we extract the form factors at the different lattice spacings and quark masses from 
nonperturbative Euclidean correlation functions. The extrapolations of these results to the physical limit will be 
discussed in Sec. V. 

We use the following interpolating fields for the At,, A c , and the proton, 

A ba = e abc d% u b b c a , (26) 

A c « = e abc (C 75 )/j 7 d% u b c c a , (27) 

N a = e abc {C lb ) M u% dh u c a , (28) 

where C is the charge-conjugation matrix, a, b , c are color indices, and a, /?, 7 are spinor indices (the symbol N is used 
for the proton to avoid confusion with the Af,-momentum p ). The tilde on the quark fields indicates gauge-covariant 
Gaussian smearing. For the u and d quarks, the smearing parameters are the same as in Ref. [61]. In the notation of 
Ref. [61], for the charm quarks we used (a, ns) = (3.0,70) at the coarse lattice spacing and (cr, ns) = (4.0,70) at the 
fine lattice spacing, and for the bottom quarks ( a , ns) = (2.0,10) at the coarse lattice spacing and (cr, ns) = (2.67,10) 
at the fine lattice spacing. The smearing of both the charm and bottom quark fields was done using Stout-smeared 
gauge links [62] with ten iterations and staple weight p = 0.08 in the spatial directions. 

In the following, we denote the final-state interpolating field by X a (= N a , A ca ) and the renormalized currents as 
Jr, where 


^7 A — Ko 
T —A 


(29) 

(30) 
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FIG. 2. Illustration of the quark field contractions on a given background gauge field for the forward (left) and backward (right) 
three-point functions. The u, d, and c quark propagators are common to the forward and backward three-point functions and 
have a Gaussian-smeared source at (xo,x). We sum over the spatial points x and y with the appropriate phases to project to 
definite momenta. The b-quark propagators are computed using the sequential source method, with sequential sources on the 
time slices xo ± t. 


with V : , and A M given by Eqs. (18-21). We set the A& three-momentum p to zero, and compute “forward” and 


“backward” three-point functions (where t > if > 0), 

c '«a’ fw) ( r > p'iMO = ^2 e ~ ip ix ~ y) (x s {x 0 ,x.) Jl(x 0 -t + t',y) A ba (x 0 - t,z)^, (31) 

C^\T, =E e ' iP ' (y_X) ( A ^( a:o +^ z ) Jr{x 0 + t',y) X s (x 0 ,x)y (32) 

y z 

as well as the two-point functions 

C 8a X,tw) (v',t) = e ~ ip ' iy ~ x) (Xs(x 0 + t,y) X a (x 0 ,x)) , (33) 

y 

Csa X ’ bv,) (p',t) = ^e _ip,(x ” y) (X,s(a;o,x) X a (x 0 -t, y)) , (34) 

y 

c Sa Ab ’ U ' ) ( t ) =^ 2 (A b s{xo + t,y) A bQ (a;o,x)) , ( 35 ) 

y 

c (2,A b ,bw)^ _ ^Af, 5 (xo, x) A ba (xo — t, y)). (36) 

y 


These definitions are similar to those in the static 6-quark case [29, 53], but with the relativistic heavy-quark action 
used here, the 6 quark can propagate in all directions, and we included additional sums over the spatial coordinates for 
the momentum projections. The quark-field contractions for the three-point functions are illustrated in Fig. 2. Only 
the 6-quark sequential propagators need to be recomputed for each source-sink separation, t. For the proton final 
state, 16 times as many sequential propagators are needed as for the A c final state because of the different structure 
of diquark contractions. The 6-quark propagators decay extremely fast with distance, and care has to be taken to 
perform sufficiently many conjugate-gradient iterations to get an accurate solution up to the distance needed. 

We computed the three-point functions for all final-state momenta p' with |p'| 2 < 12 (27r/L) 2 , and for the ranges of 
source-sink separations shown in Table V. In a first run we computed the three-point functions for all possible values 
of t/a in the wide ranges shown in the left column of Table V, but only for the lattice currents of the form qTb and 
qT 7 j j 6. In a second run, we then computed the three-point functions for all of the remaining (H(a)-improvement 
currents shown in Eqs. (18-21), but only for the subsets of separations in the right column of Table V to save computer 
time and disk space. For one of the data sets (C14), we performed the calculation of all the currents for the whole 
range of source-sink separations. As shown in Fig. 4, the effects of the additional 0(a) improvements are small. Our 
method for effectively including these corrections for all source-sink separations will be explained further below. 

To discuss the spectral decomposition of the correlation functions, we introduce the overlap factors 

(0|A 6a (0)|A 6 (p, s)) = [(Z™ + 4 2 fc V) «(p, s)] Q , (37) 

<0|X Q (0)|X(f/, S ')> = [(4 1} + 4V) U(p',s')\ a . (38) 

The two separate Z factors for each matrix element are needed because the spatial-only smearing of the quark fields 
in the interpolating field breaks hypercubic symmetry [33]. Because we set p = 0, we can write 


(0|A bQ (0)|A p (p, s)) = Z Ab u(p,s) a , 


(39) 










Set Partial (P(a)-improvement Full (b(a)-improvement 


C14 

t/a 

= 4.. 

15 

t/a = 4 ... 15 

C24 

t/a 

= 4.. 

15 

t/a = 5,8,11 

C54 

t/a 

= 4.. 

15 

t/a = 5,8,11 

F23 

t/a 

= 5.. 

15 

t/a = 6,10,14 

F43 

t/a 

= 5.. 

15 

t/a = 6,10,14 

F63 

t/a 

= 5.. 

17 

t/a = 6,10,14 


TABLE V. Source-sink separations used for the three-point functions for each data set. For the separations in the column “full 
(b(a)-improvement”, we computed the three-point functions for all of the (D(a)-corrections in Eqs. (18-21). For the separations 
in the column “partial (U(a)-improvement”, we computed only the corrections with coefficients Cy and c A ^. As explained 
in the main text and illustrated in Fig. 4, the effects of the missing terms are very small and practically independent of the 
source-sink separation, and we achieve full (b(a)-improvement for all separations by applying t-independent correction factors 
computed using the subsets of separations were all (b(a)-corrections are available. 


where Z Ab = + Z x \ Further, we introduce the following short-hand notation for the form factor decomposition 

of the matrix elements (cf. Sec. II): 

(X(p',s')\J r \A b (p,s)) = u x (p',s') &[T] u Ab (p,s) . (40) 

The spectral decompositions of the correlation functions then read 

(4' } + z x\ 0 )( m x + /) ^[r] m Ab {l +7°) 

(41) 

m Ab ( 1 + 7 °) W\ (mx +/)( 4 1} + 4 2 ) 7 °) 

(42) 


C ( 3 ’ fw) (p', i\m')= z A 

+(e 

C ( 3 ’ bw V, T,t,t-t') = Z A 


1 1 


e -Ex(t-t') e ~m Ab t' 


2E X 2m Ab 
+ (excited-state contributions), 

^ ^ -Ext' 

b 2 E x 2m Ab 
+ (excited-state contributions), 


C (2 ,- Y ’fw)( p /,f) = C^ x ’ hw \p\t) = 


1 


0 —Ext 


2 E x 

+ (excited-state contributions), 

C (2,A i ,,fw)^ = C-(2,A b ,bw)^ = _*_ e -m Ab t + 7°)] 

b 

+ (excited-state contributions), 


( 4 } + 4 2) 7 °)(m x +/)( 4 ) + 4 2 ) 7 °) 


(43) 

(44) 


where S?[r] = 7o^[r]^7oj and all correlators are 4x4 matrices in spinor space. In the above expressions, we have 
explicitly shown only the ground-state contributions, which correspond to the positive-parity baryons of interest. The 
excited-state contributions decay exponentially faster with the time separations t' , f. 

To extract individual form factors, we contract the currents in the three-point functions with suitable polarization 
vectors and form certain double ratios that eliminate all the time-dependence and overlap factors for the ground-state 
contributions. For an arbitrary four-vector n , we define 


r[n ] = n — 


0 q-n ) 


(45) 


where q = p — p' is the four-momentum transfer. By construction, r[n] is orthogonal to q. For the vector current, we 
define the three ratios 




r M [(l,0)] r„[(l,0)]Tr 

C( 3 ’ fw )(p', Y, t,f) c^’ hw \ p', 7 v ,t,t-t’) 

Tr 

C ( 2,X,av)(p/ jt ) 

Tr 

C(2,A b ,av)( f ) 



r M [(0, e j x p')] rv[(0, efc x p')] Tr 

C r(3,fw)( p / j 7 M,t,t')757 J C (3 -bw)(p', r ) u , t, t — t')757 fc 

Tr 

C'(2, X ,av)(p' ; t) 

Tr 




(46) 

(47) 


q,fj, q» Tr 

C'( 3 ’ fw )( P / , c( 3 ’ bw i( P ', y, 

t,t — t') 

Tr 

C( 2,X,av)(p^ )t ) 

Tr 




•4 (p'.mO 


(48) 
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Set 

am Ab 

am Aa 

am n 

ams c 

ams 

C14 

3.305(11) 

1.3499(51) 

0.6184(76) 

3.60327(42) 

3.0649(27) 

C24 

3.299(10) 

1.3526(57) 

0.6259(57) 

3.60312(45) 

3.0628(29) 

C54 

3.3161(71) 

1.3706(40) 

0.6580(39) 

3.60326(44) 

3.0638(33) 

F23 

2.469(16) 

1.008(12) 

0.4510(86) 

2.73156(44) 

2.3198(32) 

F43 

2.492(11) 

1.0185(67) 

0.4705(42) 

2.73169(44) 

2.3230(26) 

F63 

2.5089(70) 

1.0314(40) 

0.5004(25) 

2.73257(33) 

2.3221(22) 


TABLE VI. Hadron masses in lattice units. 


where ej is the three-dimensional unit vector in j-direction, and x is the three-dimensional vector cross product. We 

sum over repeated indices //, v from 0 to 3 and over repeated indices j,k from 1 to 3. The quantities C < - 2,Ab,av ' > and 
C (2,.Y,av) 

in the denominators are the averages of the forward- and backward two-point functions. 

These ratios are designed to isolate particular helicity form factors and are equal to 


^+(P \M') = 


(E x - m x ) 2 (E x + m x ) m Ab (m Ab + rn x )f+ 


1 2 


&±(p’,t,t) = 


4m\ b E x q 4 
(E x - m x ) 2 (E x + m x ) m Ab f± 


+ (excited-state contributions), 


* V (p',t,t') = 


m A b E x 

(E x + m x ) m Ab (m Ab -m x )fo 


+ (excited-state contributions), 


1 2 


4 Exmi 


+ (excited-state contributions). 


(49) 

(50) 

(51) 


Sample numerical results for and are shown in Fig. 3. We further define the quantities 


Rf+( Ip'I.*) 

RfoilP'lt) 


2 q 2 


(E x - m x )(m Ab + m x ) V E x + m x 


Ex 


-&+(\p'\,t,t/2), 


1 


E x 


E x - m x V E x + m x 


Ex 


m Ab - m x V E x + m x 


:^I(|p'l ,t,t/2), 


(52) 

(53) 

(54) 


where we evaluate the ratios & at t' = t /2 to minimize excited-state contamination at a given value of the source-sink 
separation t (if t/a is odd, we average over i! = (t + a )/2 and t' = (t — a )/2 instead). The notation with the 
absolute value indicates that we average over the directions of p'. Equations (52-54) yield 


Rf + (\p'\,t) = / + + (excited-state contributions), (55) 

i? /x (|p , |,t) = /_l + (excited-state contributions), (56) 

Rf 0 (\p'\,t) = fo + (excited-state contributions), (57) 


where the excited-state contributions decay exponentially with t. We checked that the helicity form factors (plus the 
corresponding excited-state contributions, for the separations we utilize) are all positive by analyzing individual three- 
point functions, so that the square roots in Eqs. (52-54) give the correct signs. Although not explicitly annotated, 
the form factors in all of the above expressions depend on |p'| and on the lattice parameters. For the axial-vector 
current, we define ± 0 as in Eq. (46-48) but with T = 7^75 in the three-point functions. The axial-vector helicity 
form factors are then extracted as 


R g+ (\p'\,t) 

Rg±(.\p'\,t) 


2 q 2 


(E x + mj)K - m x ) V E x - m x 


E x 


^+(|p'|,M/ 2 ), 


E x 


E x + m x V E x - m x ~ _L 




m Ab + m x V E x - m x 


-&oi\p'\,t,t/2). 


(58) 

(59) 


RgodP'l’t) 


(60) 
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t/a = 6 




(r; H 

CS 

<N 

a 


0.006 

0.005 

0.004 

0.003 

0.002 

0.001 

0.10 

0.08 

0.06 

0.04 

0.02 

3.0 

2.5 

2.0 


1.5 r 
1.0 - 
0.5 - 


0.025 p - 
0.020 - 
0.015 - 
0.010 - 
0.005 - 


"JH 


a 

I 


0.30 
0.25 
0.20 
0.15 - 
0.10 - 
0.05 


1.01 - 

0.8 L 
0 . 6 - . 
0.4 - 
0.2 - » 


J Ai —i A c J Aj —>■ p 


t/a = 8 


t/a = 10 



0.0 


0.2 0.4 

t' (fm) 


_L_ 

0.6 0.0 0.2 


I i i i i i 



0.4 0.6 

t' (fm) 


0.8 0.0 0.2 


0.4 0.6 0.8 

t! (fm) 


1.0 


FIG. 3. Numerical results for the vector-current ratios (46), (47), (48) and their axial-vector counterparts, at p'| 2 = 3(27r/L) 2 , 
plotted for three different source-sink separations t. The data shown here are from the C24 data set. 


When evaluating the ratios, we take the baryon masses in lattice units, am \ b , arriA c , and amjv, from exponential fits 
to the zero-momentum two-point functions for each data set; see Table VI. We then compute the energies a£?A c (p'), 
and clEn(p') from these masses using the relativistic continuum dispersion relation, and we also compute a 2 q 2 from 
these masses and energies. Because the form factors are dimensionless, the values of the lattice spacing are not needed 
at this stage. The ratios are evaluated using statistical bootstrap, and we use corresponding bootstrap samples for 
the masses to take into account all correlations. 

As mentioned earlier, except in the case of the C14 data set, we have “full-O(a) improvement” (“FI”) data only 
for three source-sink separations in each data set, but we have data with “partial C ) (a)-improvement” (“PI”) for all 
source-sink separations in the ranges shown in Table V. To account for this, we computed the ratios 


< n (|p'l,*) 

#/ PI) (Ip'I^)’ 


(61) 


7lC'T'\ 

where / = /+, f±, /o, g+, g±, go, for those source-sink separations where both Rj ' and Rj are available. Numerical 
results for Eq. (61) from the C14 data set (where we have FI data for all values of t) are shown in Fig. 4. In the case 
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FIG. 4. Ratios of the “fully 0(a)-improved” (“FI”) and “partially 0(a)-improved” (“PI”) data for Rf{t) for the six different 
helicity form factors, at |p'| 2 = 3(2tt/Z,) 2 , from the C14 data set. The partially improved data only include the currents qT^j^jb 
for the 0(a)-improvement, where q = u,c. For A;, —> p , these are the only currents needed at tree level and the ratio is very 
close to 1. For A b — > A c , the ratio deviates from 1 significantly more, because the currents j^jTb are missing the partially 
improved data, but are needed already at tree level in this case. The range of source-sink separations shown for Af, — > p is 
smaller because the statistical fluctuations in the correlators were too large to reliably compute the individual quantities Rf(t) 
for t > 1.0 fm in this case. 


of At —► A c (Fig. 4 right), the correction is as large as 2% for some of the form factors, but is independent of the 
source-sink separation to a high degree (even though Rj FI) and Rf P1> individually have a strong t-dependence). The 
same behavior is found at other values of the momentum. In the case of Ab —> p (Fig. 4 left), the correction shows a 
more significant dependence on the source-sink separation, but is smaller than 0.3% for all form factors. For the data 
sets other than C14 we therefore performed constant fits to the ratios (61) as a function of t , individually for each 
form factor /, each momentum |p'|, and each data set. If these fits had a poor x 2 /dof, we excluded the shortest or 
the two shortest separations. In this way, we obtained correction factors, which we then applied to Rj ; (|p'|, t) at all 
separations, to effectively obtain R'j P (|p^|, t) at all separations. This procedure is accurate to better than permille 
level. In the following, all ratios i?/(|p'|,f) are understood to be corrected using this procedure for all source-sink 
separations. 

For the further data analysis, we then also formed the linear combinations 


R t? 

Rfv 
J 2 

Rfv 

R St 

R f£ 


(m Ab + m x ) 2 Rf + - q 2 Rf ± 

5 

m Ab (m Ab + m x )(Rf ± - Rf + ) 

1 

s+ 

m Ab (m Ab - m x ) [(m Ai , + m x ) 2 {Rf 0 - Rf + ) + q 2 (Rf ± - #/ 0 )] 

q 2 s+ 

[m Ab - m x ) 2 R g+ - q 2 R g _l 

? 

S_ 

m Ab {m Ab - m x ){R g+ - R g J 
s_ 

m Ab (m Ab + m x ) [(m At - m x ) 2 (R g+ - R go ) + q 2 (R go - R g± )] 

q 2 S- 


(62) 

(63) 

(64) 

(65) 

( 66 ) 
(67) 


which, according to the relations in Eqs. (8)-(13), become equal to the Weinberg form factors at large t. To extract the 
ground-state form factors from Rf(\p'\,t) (for both the helicity and Weinberg form factors), we performed correlated 
fits of the f-dependence including exponential correction terms to account for the leading excited-state contributions, 
thereby extrapolating Rf(\p'\, t) to t = oo. To discuss these extrapolations in more detail, it is convenient to denote 
the data for Rf(\p'\,t) by 


Rf,i,n(t ), (68) 

where / = /+, f±, f 0 , g+, g±, go, f¥,fX, /(S / 2 , f£ labels the form factors, i = C14, C24, C54, F23, F43, F63 labels 
the data set (cf. Table I), and n labels the final-state momentum via |p'| 2 = n (2n) 2 /L 2 . We performed the fits using 
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I A5 —> A c } A & —> p 



t (fm) 


t (fm) 


FIG. 5. Extrapolations of _R/(|p'|,f) to infinite source-sink separation. The data shown here are at momentum |p'| 2 = 
3(27r/L) 2 , and are from the C24 data set. For each momentum, all vector (or axial vector) form factors from all data sets are 
fitted simultaneously as explained in the main text. 
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the functions 


Rf,i,n(f) fi,n V ^ 


-5, 




= s n 


GeV, 


(69) 


with parameters Af^ tn , and where /q n are the form factors we aim to extract. By writing the energy-gaps 

3f,i,n in the above form, we impose the constraint <5/y, n > <5 m in- We chose <5 m j n = 170 MeV, which is smaller than 
any expected energy gap (given our prior knowledge of the hadron spectrum at our values of the pion masses). This 
constraint has negligible effect in most cases, but prevents numerical instabilities for some form factors at certain 
momenta where the data show no discernible excited-state contamination. 

At each momentum n, we perform one coupled fit to all the data for the vector-current form factors (/+, f±, /o, 
fY , fY, fY ) an d another coupled fit to all the data for the axial-vector-current form factors ( g + , g±, g 0 , fY, fY , fY)- 
This allows us to implement two additional constraints to stabilize the fits, based on the following knowledge: 

• Because the lattice size, L (in physical units), is equal within uncertainties for all data sets ( L « 2.7 fm), the 
squared momentum |p'| 2 = n{2n/L) 2 for a given n is also equal within uncertainties for all data sets. This 
means that the energy levels, and hence the parameters lf,i, n , are expected to be approximately equal across all 
data sets i, up to some dependence on the pion mass and the lattice spacing. 

• By construction, the data Rf in (t) for the helicity and Weinberg form factors exactly satisfy the defining relations 
(62)-(62) at each value of the source-sink separation. The extracted ground-state form factors fi, n should also 
satisfy these relations. 


For the coupled fit to all vector form factor data at a given momentum n, we therefore add the following terms, 
corresponding to Gaussian priors, to the x 2 function: 


Xv,r 


Xv,n+X! 

/ 


(^/,C14,n — lf,C24,n) 2 , (lf,C24,n — If, C54,n ) 2 ^ ( lf,¥23,n ~ ^/,F43,ra ) 2 




C14,C24i 2 


CT, 


C24,C54i 2 


\o, 


F23,F43-|2 


G/,F43 ,n — lf,F63,n) 2 ( lf,C54,n — lf,^3,n) 2 ^ 


cr, 


F43,F63i 2 


r C54,F63'|2 i 9 
\Om r + oi 


+ E ( /+.«." - fli,n - - -- -jLn ) /'<*) 

\ rnA b ,i(mA b ,i + m x ,i) / J 


A b,i 


E f n _ rV _ m A b ,i + mx,i f y 

I /_L,z,n Jl,i t n 
i ^ 

E! ( fo,i,n ~ fY,i,n ~ 


1°) 


(fi 

Vi.n 


-7- sf¥,i,n 

m Ab ,i{rnA b ,i - m x ,i) ” / J 


(70) 


where cr a = 0.1 and 


KE=^[(<) 2 -(mi) 2 ] 2 , (71) 

with w m = 4GeV -2 . With these widths, the first two lines in Eq. (70) implement the constraint that the energy gaps 
(<5/,i, n — £min) at given momentum n should not change by more than 10% when going from the fine to the coarse 
lattice spacing and not more than 400% times the change in m 2 (in GeV 2 ); both are reasonable assumptions given 
the prior experience with hadron spectroscopy in lattice QCD. Note that absolute variations of lf.i, n translate to 
relative variations of — <5min) because d[exp(Z/ j j iTl )]/exp (lf,i, n ) = d lf,i, n - The last three lines in Eq. (70) enforce 
the relations (62)-(64) between the ground-state vector form factors in the helicity and Weinberg definitions (we set 
07 = 10 -4 ). For the fit to the axial vector form factor data, analogous terms are added to Xa n- 

We initially included all available values of t in the fits, and then removed data points for each form factor at 
the smallest t until the fits had good quality as determined by the correlated x 2 /dof. To estimate the remaining 
systematic uncertainties associated with higher excited states, we then further removed the next-lowest values of t 
simultaneously for all Rf,i, n and computed the resulting shifts in /, „. We then took the larger of the following two 
as our estimate of the excited-state systematic uncertainty: i) the shift in f ln at the given momentum n, and ii) 
the average of the shifts fi^ n over all momenta n. We added these excited-state uncertainties in quadrature to the 
statistical uncertainties in fi, n . All result for /).„ are listed in Appendix A. As can be seen in Tables XIV and XVI, 
the results for the second-class form factor fY are very close to or consistent with zero for both A& —t A c and At, —> A c , 
despite the rather large mass differences mb — m u and Wf, — m c . The results for the other second-class form factor 
fY are significantly nonzero, but are still noticeably smaller than the results for the first-class form factors. 
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V. CHIRAL/CONTINUUM/KINEMATIC EXTRAPOLATION OF THE FORM FACTORS 

The last step in the data analysis is to perform fits of form factor results using suitable functions describing 

the dependence on the momentum transfer, the dependence on the up and down quark masses (or equivalently the 
pion mass), and the dependence on the lattice spacing. We perform global fits of the helicity form factors based on 
the simplified ^-expansion [37], modified to account for pion-mass and lattice-spacing dependence. The expansion 
parameter z* for a form factor / is defined as 



where we choose 


(72) 


to = (m At - m. x) 2 , 


(73) 


so that the point z = 0 corresponds to q 2 = q^ ax (i.e. p' = 0 in the A b rest frame). The values of t/ + are discussed 
further below. After factoring out the leading pole contribution, we expand the form factors in a power series in z*. 
We find that our lattice data can be described well by keeping only the zeroth and first order in z *. As explained 
further below, we also perform higher-order fits to estimate systematic uncertainties. Our nominal (as opposed to 
higher-order) fits are of the form 


/(<? 2 ) = 


1 


1 - q 2 /(m f p oiJ 2 
,'|2 


X 


1 + b f 


(tt/o ) 2 


J '0 l 1 + c o 


(] f A QCd 


j m l -m 2 phys 


Al 


+ a{ z f (q 2 ) 


(tt/u) 2 ]’ 


(74) 


with fit parameters <Zg, a{, Cq, b?, and dA Here, m v are the valence pion masses of each data set (see Table I), and 
w-n-.phys = 134.8 MeV is the physical pion mass in the isospin limit [63]. As discussed in Ref. [29], chiral-perturbation- 
theory predictions for the pion-mass dependence of the form factors considered here are unavailable and would be 
of limited use because of the large momentum scales in these matrix elements, and because of the large number of 
low-energy constants. In Eq. (74) we describe the pion-mass dependence through the factor 1 + Cg (m 2 — m 2 phys )/A 2 
multiplying do- Here, we introduced the scale A x = Anf with / = 132MeV so that Co becomes dimensionless. Because 
our lattice actions and currents are (!l(a)-improved, we allow for a quadratic dependence on the lattice spacing via 
the factor in the second line of Eq. (74), where Aqcd = 0.5 GeV. The parameters b? and df describe the momentum- 
dependent and momentum-independent parts of the lattice discretization errors. We use the individual lattice QCD 
results for the baryon masses from each dataset (see Table VI) to evaluate a 2 q 2 and z, and we take into account the 
uncertainties and correlations of these masses. We set the pole masses equal to 


aTO poie = amps + aA f , (75) 

where amps is the pseudoscalar B u or B c mass (in lattice units) computed individually for each data set (and also 
listed in Table VI), and A* is the mass splitting (in GeV) between the meson with the relevant quantum numbers 
and the pseudoscalar B u (for A b —► p) or B c (for A b —► A c ). We use fixed values of for all data sets, based on 
experimental data (where available) [1] and averages of our lattice QCD results over the different data sets. These 
values are given in Table VII. The pole factor is then written as 


1 — (a 2 q 2 )/(amps + aAf)' 2 ’ 

so that the explicit value of the lattice spacing is needed only for the term a AC Note that when the input values of 
A / are varied, the shape parameters Oq and a{ returned from the fit change in such a way as to largely cancel the 
effect of this variation on the form factors (varying A* by 10% changes the form factors themselves by less than 1%). 

f 

The parameter t\ should be set equal to or below the location of any singularities remaining after factoring out 
the leading pole contribution. The A b — > p and A b -A A c form factors (in infinite volume) have branch cuts starting 
at q 2 = (rriB + rn w ) 2 and q 2 = (m^ + m^) 2 , respectively. In the case of A^ —» p, the form factors /o, g+, g± have no 
poles below this branch cut, and the form factors /+, f±, go only have a single pole below q 2 = (ms + m w ) 2 , which 
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/ 

j p 

t+( Ai> -> P ) 

TOpol e (A6 ->■ P) 

A f (A b ->• p) 

—>• A c ) 

TO-pole (A& ->• A c ) 

A f (A b —>■ A c ) 

f+, f± 

I - 

(rriB + ran) 2 

mB + A* 

46 MeV 

( m pole) 2 

mB c + A f 

56 MeV 

fo 

0+ 

(mB + rn n ) 2 

mB + 

377 MeV 

( m pole) 2 

niB c + A f 

449 MeV 

g+, g± 

i+ 

(mB +m n ) 2 

mB + A* 

427 MeV 

( m pole) 2 

mB c + A f 

492 MeV 

go 

(T 

(mB + nin) 2 

mB + A* 

0 

( m pole) 2 

n%B a + A f 

0 


TABLE VII. Values of v + and m^ ole . To evaluate the form factors in the physical limit, ms = 5.279 GeV, m n = 134.8 MeV, 
and ms c = 6.276 GeV should be used. 


gets removed by the explicit factor of 1/[1 — <? 2 /(trip ole ) 2 ]- We therefore set t+ = (tob + to w ) 2 for all A b —> p form 
factors. More precisely, to evaluate the dimensionless quantity a 2 tl in the fit, we use 


a 2 t f + = (am PS + arn^,phys) 2 (for A b -S > p), (77) 

where amps is the pseudoscalar B u mass (in lattice units) computed individually for each data set, and m^phys = 
134.8 MeV. This means that the value of the lattice spacing is needed only for the term a?n w , p h ys - We checked that 
using the individual lattice pion masses of each data set instead of the physical pion mass has a negligible effect on 
the extrapolated form factors. In the case of A b —> A c , the onset of the branch cut, tob + Tod, is several hundred 
MeV above the lowest pole for all form factors, and there may be additional poles below ms + m,D- We therefore set 
t+ = (TOp oIe ) 2 for the A b —» A c form factors; more precisely, 

a 2 t+ = (amps + aA^) 2 (for A b —► A c ), (78) 


where amps is the pseudoscalar B c mass (in lattice units) computed individually for each data set, as discussed above. 
With this choice of t+, the factors of 1/[1 — g 2 /(m.p ole ) 2 ] are not strictly necessary, but we find that they improve the 
quality of the fit at first order in the z expansion. 

We implement the constraint g±(q^ nax ) = g+(t 7max) [Ed- (16)] at z 9±,9+ = 0 and a = 0 by using shared parameters 
al ±:9+ and Cg~ L,5+ for these two form factors. We impose the constraints /o(0) = /+(0) and go(0) = g + ( 0) [Eqs. (14) 
and (15)] using Gaussian priors with widths equal to max^ 0 (0), z i+ (0)] 2 and ma,x[z 90 (0), z 9+ (0)] 2 , respectively, to 
allow for missing higher-order terms in z*. For A b —> p , we performed one global fit to all helicity form factors, taking 
into account the correlations between different form factors, different momenta, and different data sets. For A b —> A c , 
such a global fit showed indications of problems associated with a poorly conditioned data covariance matrix, and 
we additionally performed fits of the subsets {/+,/o}, {f±}, {ff+:ff_L,5o} to reduce the sizes of the data covariance 
matrices. We then took the central values and covariances of the form factor parameters within each subset from 
these subset fits, and only used the global fit to estimate the cross-covariances between the parameters in different 
subsets. 

The physical limit is given by a —> 0 and —» m^phys, and correspondingly Eq. (74) reduces to the simple form 


f(4) 


1 — q 2 /(m: 


pole 


[4 + 4 *V)] 


(79) 


where q 2 should be evaluated using the experimental values of the baryon masses, and m^ ole , t+ should be set to 
the values given in Table VII, with tob = 5.279 GeV, to n = 134.8 MeV, tob c = 6.276 GeV. The central values and 
uncertainties of the parameters {ag, a{} from the nominal fit are given in Table VIII, and the correlation matrices 
are given in Table IX. The parameter covariances cov(p, q) can be obtained from the correlations corr(p, q) and 
uncertainties a p , a q using cav(p , q) = a p a q corr(p, g); the central values and covariance matrices of the fit parameters 
are also provided as ancillary files with the arXiv submission of this article. Plots of the lattice data along with the 
physical-limit fit curves are shown in Figs. 6, 7, 8, and 9. 

To estimate the systematic uncertainties caused by our assumptions on the lattice-spacing, quark-mass, and q 2 - 
dependence, we also perform fits that include additional higher-order terms, employing the form 


/ho(<7 2 


1 


1 - g 2 /(TOp 0le ) 2 0 

+a{ [^(g 2 )] 2 l \l + b f 


f ( / to 2 -to 2 h _ f 

-0 1 + 4 - A2 + 'o 


m; — m 


^phys 


V12 


(n/a) 2 


d f a qcd 


b f 


-*./13 


(-7T/a) 2 ( n/a ) 3 


J/ A QCD + 
(n /a) 3 


( 1 + 4 m '- A f*- )^tf) 

IP'IAJ 


■f Ip'I 2a qcd 
(n/a) 3 


- k f - 


QCD 


(tt /a) 3 


■ (80) 
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Parameter 

A b ->• p 

A b —y A c 

a f 0 + 

0.4382 ± 0.0315 

0.8146 ±0.0167 

a{+ 

-0.6452 ± 0.2093 

-4.8990 ± 0.5425 

a f 0 ° 

0.4189 ± 0.0256 

0.7439 ±0.0125 

a{° 

-0.7862 ±0.2038 

-4.6480 ± 0.6084 

a o ± 

0.5389 ± 0.0435 

1.0780 ± 0.0256 

a f y 

-0.8069 ± 0.3039 

-6.4170 ± 0.8480 

a.r 9+ 

0.3912 ±0.0198 

0.6847 ±0.0086 


-0.8167 ±0.1749 

-4.4310 ±0.3572 

a 90 

a 0 

0.4526 ± 0.0292 

0.7396 ± 0.0143 

af 

-0.7817 ±0.1886 

-4.3660 ± 0.3314 

a i ± 

-0.9061 ±0.1956 

-4.4630 ± 0.3613 


TABLE VIII. Central values and uncertainties of the nominal form factor parameters for At —> p and A;, —> A c . See Table IX 
for the correlation matrices. 


This allows for higher-order variation in the lattice spacing, quark masses, and momentum dependence. The data 
themselves do not determine this more complex form sufficiently well, so we constrain the higher-order coefficients 
Cq, c{, 6f d,f , jf, k,f to be natural-sized using Gaussian priors with central value 0 and width 10. We constrain the 
second-order z-expansion coefficients a J 2 using Gaussian priors with central values 0 and widths given by approximately 
twice the magnitude of the previous (nominal) fit results for a{. Given that this fit is quadratic in z? . we now impose 
the kinematic constraints (14) and (15) at q 2 = 0 up to widths of max[z/°(0), z?+ (0)] 3 and max[z 90 (0), z 9+ (0)] 3 , 
respectively. 

In the higher-order fit, we use bootstrap data for the correlator ratios in which the matching- and 0(a)-improvement 
coefficients were drawn from Gaussian random distributions with central values and widths according to Table III. 
Thus, the higher-order fit results also include the perturbation-theory systematic uncertainty. To take into account 
the uncertainties of the lattice spacings, we promote the lattice spacings of the different ensembles to fit parameters, 
constrained with Gaussian priors according to the central values and uncertainties given in Table I. The systematic 
uncertainties caused by the finite lattice volume cannot easily be estimated from the data, because all of our data 
sets have approximately the same lattice size, L ss 2.7 fm. Finite-volume effects have been calculated using chiral 
perturbation theory for the nucleon magnetic moment [64] and axial charge [65], and, specifically for the ensembles 
used herein, for the heavy-baryon axial couplings [61, 66]. Based on this experience, we estimate that the finite-volume 
systematic uncertainties in our results are 3% for the A & —>■ p form factors, and 1.5% for the A^ —> A c form factors. 
The neglected isospin breaking effects in the form factors are estimated to be of order 0((rrid — to m )/Aqcd) ~ 0.5% 
and 0(a e , m ,) ~ 0.7%. Finally, there is an uncertainty resulting from the tuning of the relativistic heavy quark (RHQ) 
parameters, which was performed in Ref. [51] for the 6 quark and in Ref. [52] for the c quark. In Ref. [12], the same 
6-quark parameters were used to compute the B —> ir form factors on the same gauge field configurations as in the 
present work, and the uncertainties of am^ , and c^ BB were propagated to the form factors by repeating the 
calculation for multiple values of awg\ ^ h> , and c B ^ B . The resulting uncertainties in the B — > n form factors were 
found to be 1%. We could not afford to repeat the present calculation for multiple values of the RHQ parameters, 
and therefore adopt the 1% estimate also for the 6-quark parameter uncertainty in the A& —> p and A;, —> A c form 
factors. We are unable to estimate the c-quark parameter uncertainty in the A b —> A c form factors at this time, but 
we note that our choice of parameters precisely reproduces the experimental values of the charmonium masses and 
hyperfine splittings, as well as the relativistic continuum dispersion relation, on both the coarse and the fine lattices 
[52]. To estimate the effect of the light-quark-mass uncertainties, we promote the pion masses (in lattice units) to 
fit parameters, constrained with Gaussian priors according to the central values and widths given in Ref. [36]. We 
find that this results in a smaller than 0.1% uncertainty in the A b —> p and A b —> A c form factors. To incorporate 
the finite-volume, isospin breaking, and RHQ parameter tuning uncertainties in the higher-order fit, we add these 
uncertainties to the data covariance matrix before performing the fit. We assume that these uncertainties are 100% 
correlated between the different data sets and different final-state momenta, and between the three different form 
factors corresponding to the same type of current (vector or axial vector). 

In the physical limit, the higher-order fit functions reduce to 

/ ho (<? 2 ) = - - ) f [a f 0 + a{ z f (q 2 ) + a{ [z f (q 2 )} 2 ] . (81) 

1 - q 2 /{m J pole y 
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The values of the parameters Oq, a{, a^, their total uncertainties, and their correlation matrices are given in Tables X 
and XI, and are also included as ancillary files with the arXiv submission. The recommended procedure for computing 
the central value, statistical uncertainty, and total systematic uncertainty of a general observable depending on the 
form factor parameters (for example, a differential decay rate at a particular value of q 2 , or an integrated decay rate, 
or a ratio of decay rates) is the following: 

1. Compute the observable and its uncertainty using the nominal form factors given by Eq. (79), with the parameter 
values and correlation matrices from Tables VIII and IX. Denote the so-obtained central value and uncertainty 
as 

O, a 0 . (82) 

2. Compute the same observable and its uncertainty using the higher-order form factors given by Eq. (81), with 
the parameter values and correlation matrices from Tables X and XI. Denote the so-obtained central value and 
uncertainty as 

Oho, co,ho- (83) 


3. The final result for the observable is then given by 


O ± <7 0 ± max (|Oho - O |, yJWo^o ~ a o\) • 


syst. 


(84) 


In other words, the central value and statistical uncertainty are obtained from the nominal fit, and the systematic 
uncertainty is given by the larger of the following two quantities: i) the shift in the central value between the nominal 
fit and the higher-order fit, and ii) the increase in the uncertainty (computed in quadrature as shown above) from 
the nominal fit to the higher-order fit. The statistical and systematic uncertainties in Eq. (84) should be added in 
quadrature. By construction, the above procedure gives the combined systematic uncertainty associated with the 
continuum extrapolation, chiral extrapolation, z expansion, perturbative matching, scale setting, 6-quark parameter 
tuning, finite volume, and missing isospin symmetry breaking/QED. 

Plots of the form factors including the systematic uncertainties, computed as explained above, are shown in Figs. 10 
and 12. The relative systematic uncertainties in the form factors are shown in Figs. 11 and 13. In addition to the 
combined systematic uncertainty (thick black curves), these figures also show the individual sources of uncertainty. 
The individual systematic uncertainties were estimated using additional fits as follows: 

• Continuum extrapolation uncertainty: only the higher-order terms with coefficients 6-Q <V, jf, kf were added 
to Eq. (74). 

• Chiral extrapolation uncertainty: only the higher-order terms with coefficients q), c{ were added to Eq. (74). 

• 2 expansion uncertainty: only the higher-order term a{ [z? (q 2 )] 2 was added to Eq. (74). 

• Matching & improvement uncertainty: no higher-order terms were added to Eq. (74), but the the matching- and 
©(^-improvement coefficients were drawn from Gaussian random distributions with central values and widths 
according to Table III when computing the correlator ratios using bootstrap. 

• Scale setting (i.e., lattice spacing) uncertainty: no higher-order terms were added to Eq. (74), but the lattice 
spacings were promoted to fit parameters constrained with Gaussian priors according to the central values and 
uncertainties given in Table I. 

• Finite-volume effects, missing isospin breaking/QED, and RHQ parameter tuning uncertainties: no higher-order 
terms were added to Eq. (74), but the estimates of these uncertainties (as discussed above) were added to the 
data covariance matrix used in the fit. 

Figures 11 and 13 show that near q 2 = <// nax , the finite-volume and chiral-extrapolation uncertainties are the largest, 
but as the momentum |p'| increases (corresponding to decreasing q 2 ), the z-expansion and continuum extrapolation 
uncertainties grow and become dominant. The continuum extrapolation uncertainty should not be interpreted as the 
actual size of lattice discretization errors; the reason for the large continuum extrapolation uncertainty is primarily 
that we have only two lattice spacings and our data do not tightly constrain all of the extrapolation coefficients. 

Discretization errors associated with the relativistic heavy quark actions used for the 6 and c quarks are not 
necessarily well approximated by the leading terms in an expansion in oAqcd and aph These errors can be described 
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by mismatches of the coefficients of higher-dimensional operators in the heavy-quark expansions of the lattice theory 
and continuum QCD [47, 67-69]. In Ref. [12], the resulting heavy-quark discretization errors in the B —► tt form 
factors were estimated using power counting to be of order 2% for the same lattice actions and parameters as used 
in the present work. For most of the kinematic range, our estimate of the total continuum-extrapolation uncertainty 
in the Ab —» p form factors is larger than this power-counting estimate, see Fig. 11. Similarly, a comparison with 
the analysis of B —» D form factors in Ref. [59] suggests that heavy-quark discretization errors in the A& —»• A c form 
factors are smaller than or compatible with our estimates of the total continuum-extrapolation uncertainties in the 
entire kinematic range. 


Aj, — > p 

4+ 

a i + 

a /° 

a 0 

a{° 

a f - l 

a o 

a i ± 

n g± ’ a + 

a Q 

a?+ 

n so 
a Q 

af 

af-L 

a o + 

1 

-0.9058 

0.5081 

-0.4403 

0.5299 

-0.3987 

0.5362 

-0.4112 

0.6302 

-0.5305 

-0.3898 

a{ + 

-0.9058 

1 

-0.4280 

0.4312 

-0.4238 

0.3739 

-0.4402 

0.3912 

-0.4901 

0.4839 

0.3668 

4° 

0.5081 

-0.4280 

1 

-0.8533 

0.4251 

-0.3226 

0.6963 

-0.5274 

0.5504 

-0.4694 

-0.4886 

a[° 

-0.4403 

0.4312 

-0.8533 

1 

-0.3525 

0.3008 

-0.5963 

0.5917 

-0.4661 

0.4554 

0.5667 

4 ± 

0.5299 

-0.4238 

0.4251 

-0.3525 

1 

-0.8930 

0.4748 

-0.3554 

0.5975 

-0.5111 

-0.3348 

a i ± 

-0.3987 

0.3739 

-0.3226 

0.3008 

-0.8930 

1 

-0.3664 

0.3156 

-0.4349 

0.4542 

0.2900 

a 0 ±,9+ 

0.5362 

-0.4402 

0.6963 

-0.5963 

0.4748 

-0.3664 

1 

-0.8434 

0.6238 

-0.5500 

-0.7905 

4 + 

-0.4112 

0.3912 

-0.5274 

0.5917 

-0.3554 

0.3156 

-0.8434 

1 

-0.4761 

0.5011 

0.8778 

n 90 

a 0 

0.6302 

-0.4901 

0.5504 

-0.4661 

0.5975 

-0.4349 

0.6238 

-0.4761 

1 

-0.9039 

-0.4497 

af 

-0.5305 

0.4839 

-0.4694 

0.4554 

-0.5111 

0.4542 

-0.5500 

0.5011 

-0.9039 

1 

0.4632 

af- 

-0.3898 

0.3668 

-0.4886 

0.5667 

-0.3348 

0.2900 

-0.7905 

0.8778 

-0.4497 

0.4632 

1 

A b — y A c 

°o + 

a i + 

O fo 
a 0 

a{° 

a 0 

4- 

n 9±,9+ 

a o 

4+ 

n 90 

a 0 

af 

af- 

ai+ 

1 

-0.6644 

0.6827 

-0.4853 

0.6218 

-0.3906 

0.4828 

-0.3152 

0.5636 

-0.4317 

-0.3763 

af+ 

-0.6644 

1 

-0.6515 

0.9445 

-0.3853 

0.5109 

-0.3831 

0.4915 

-0.2979 

0.4916 

0.4764 

4° 

0.6827 

-0.6515 

1 

-0.7040 

0.4208 

-0.3620 

0.6174 

-0.4822 

0.4320 

-0.4726 

-0.4756 

a{° 

-0.4853 

0.9445 

-0.7040 

1 

-0.2738 

0.4739 

-0.3888 

0.5261 

-0.2164 

0.4779 

0.4877 

a o x 

0.6218 

-0.3853 

0.4208 

-0.2738 

1 

-0.6637 

0.3933 

-0.2369 

0.5161 

-0.3639 

-0.2926 

a^ 

-0.3906 

0.5109 

-0.3620 

0.4739 

-0.6637 

1 

-0.2903 

0.3509 

-0.2443 

0.3640 

0.3400 

a o ± ’ 9+ 

0.4828 

-0.3831 

0.6174 

-0.3888 

0.3933 

-0.2903 

1 

-0.7304 

0.6365 

-0.6743 

-0.7301 

4 + 

-0.3152 

0.4915 

-0.4822 

0.5261 

-0.2369 

0.3509 

-0.7304 

1 

-0.3829 

0.8725 

0.9171 

a 90 

a 0 

0.5636 

-0.2979 

0.4320 

-0.2164 

0.5161 

-0.2443 

0.6365 

-0.3829 

1 

-0.6843 

-0.4846 

af 

-0.4317 

0.4916 

-0.4726 

0.4779 

-0.3639 

0.3640 

-0.6743 

0.8725 

-0.6843 

1 

0.8456 

af- 

-0.3763 

0.4764 

-0.4756 

0.4877 

-0.2926 

0.3400 

-0.7301 

0.9171 

-0.4846 

0.8456 

1 


TABLE IX. Correlation matrices of the nominal form factor parameters for A& —» p (top) and A;, —► A c (bottom). 
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h£hC14 ^C24 ^C54 ^ F23 h^F43 ^ F63 - a = 0, m v = 134.8 MeV 





FIG. 6. At, —> p vector form factors: lattice results and extrapolation to the physical limit (nominal fit). The bands indicate 
the lcr statistical uncertainty. 
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h£hC14 ^C24 ^C54 ^ F23 h^F43 ^ F63 - a = 0, m v = 134.8 MeV 





FIG. 7. Ab —> p axial-vector form factors: lattice results and extrapolation to the physical limit (nominal fit). The bands 
indicate the la statistical uncertainty. 
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|-^C14 h|hC 24 i-$h C54 ^ F23 F43 ^ F63 - a = 0, m„. = 134.8 MeV 





FIG. 8. Ab —¥ A c vector form factors: lattice results and extrapolation to the physical limit (nominal fit). The bands indicate 
the lcr statistical uncertainty. 
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h£hC14 i-J-i C24 ^C54 ^ F23 h^F 43 >-f-<F63 - a = 0, to* = 134.8 MeV 





FIG. 9. A b —> A c axial-vector form factors: lattice results and extrapolation to the physical limit (nominal fit). The bands 
indicate the lcr statistical uncertainty. 
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Parameter A& —> p At, —> A c 


a 0 + 

0.4251 

± 

0.0388 

0.8103 

± 

0.0276 

a{+ 

-0.7088 

± 

0.3361 

-4.7480 

± 

0.9429 

a f 2 + 

0.8925 

± 

0.8869 

0.7862 

± 

8.8020 

a 0° 

0.4144 

± 

0.0321 

0.7389 

± 

0.0225 

a{° 

-1.0420 

± 

0.3142 

-4.5630 

± 

0.9426 

4° 

1.9260 

± 

0.9190 

2.7050 

± 

8.4430 


0.5214 

± 

0.0520 

1.0940 

± 

0.0435 

a[- 

-0.8247 

± 

0.4424 

-6.4410 

± 

1.5010 

02 X 

0.7609 

± 

1.2770 

2.3160 

± 

11.320 

a 9 4' 9+ 

0.3889 

± 

0.0260 

0.6848 

± 

0.0184 


-1.0730 

± 

0.2617 

-4.3790 

± 

0.6954 

<4 + 

1.9860 

± 

0.8247 

1.2810 

± 

7.3650 

n 90 

a 0 

0.4419 

± 

0.0388 

0.7408 

± 

0.0258 

a{° 

-0.8649 

± 

0.3481 

-4.3860 

± 

0.8774 

a 9 2 ° 

0.9969 

± 

0.8955 

1.3380 

± 

8.0440 


-1.0840 

± 

0.2732 

-4.6270 

± 

0.7088 


1.4520 

± 

1.0680 

1.6140 

± 

7.4530 


TABLE X. Central values and uncertainties of the higher-order form factor parameters for At, —> p and At, —► A c . See Table XI 
for the correlation matrices. 







0.7671 0.2482 0.5337 -0.2670 -0.0922 0.5121 -0.2469 -0.0180 0.3774 -0.2148 -0.0472 0.4420 -0.2680 0.0018 -0.2284 -0.0231 
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TABLE XI. Correlation matrices of the higher-order form factor parameters for A b p (top) and At, —> A c (bottom). 
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q 2 (GeV 2 ) q 2 (GeV 2 ) 


FIG. 10. Final results for the Aj, —> p form factors. The inner bands show the statistical uncertainty and the outer bands show 
the total uncertainty. 
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FIG. 11. Systematic uncertainties in the Af, —> p form factors in the high-q 2 region. As explained in the main text, the combined 
uncertainty is not simply the quadratic sum of the individual uncertainties. 
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FIG. 12. Final results for the At, —> A c form factors. The inner bands show the statistical uncertainty and the outer bands 
show the total uncertainty. 
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FIG. 13. Systematic uncertainties in the A b —»■ A c form factors in the high-g 2 region. As explained in the main text, the 
combined uncertainty is not simply the quadratic sum of the individual uncertainties. 
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VI. PREDICTIONS FOR THE A b ->p£~v e AND A 6 -a A c £~u e DECAY RATES 


In this section, we present predictions for the Aft -A pl~Vi and A& -A A c l~vn differential and integrated decay rates 
using our form factor results. Including possible right-handed currents with real-valued e^, the effective Hamiltonian 
in Eq. (2) leads to the following expression for the differential decay rate in terms of the helicity form factors, 

dr _ G%\V£\ 2 y/s^sZ / m 2 \ 2 

dq 2 768ir 3 m\ b \ q 2 ) 

x 14 (ml + 2 q 2 ) (s+ [(1 - e*)g ± ] 2 + s_ [(1 ± ef )/j_] 2 ) 

+ 2 m t , +2 q ^ [( mAt - mx ) (1 - e q )g +\ 2 + s_ [{m Ab + m x ) (1 + ef)/+] 2 ) 

±^T ( s + [( mA fc “ m *) ( 1 + e q)M 2 + s - [( TO A t + to -y) (1 - ef)g 0 ] 2 ) |, (85) 

where, as before, X = p, A c denotes the final-state baryon, and 

s± = (m Ab ± m x ) 2 - I 2 - (86) 


Expressions for the individual helicity amplitudes and the angular distributions can be found in Refs. [27, 28, 70]. By 
combining experimental data with our form factor results, novel constraints in the (V^, e^) plane can be obtained. 

In the following, we consider the Standard Model with = V qb and e q = 0. Our predictions of the Aft -A p£~P( 
and Aft -A A c £~vg differential decay rates for £ = e,fi,T are shown in Figs. 14 and 15. The central values, statistical 
uncertainties, and systematic uncertainties have been calculated using Eq. (84); all baryon and lepton masses were 
taken from Ref. [1]. Our results are most precise in the high-g 2 region, where the form factor shapes are most tightly 
constrained by the lattice QCD data. We obtain the following partially integrated decay rates 


C P ^(15GeV 2 ) 


Ca c /xp( 7 GeV 2 ) 


and their ratio 


1 dr(Aft-S>p^t Vp) 2 

IWAsGeV* dg 2 q 

1 /" 3m “ dr(Aft —)• A c /r“^) 2 

\V cb \ 2 J 7Ge y 2 dq 2 Q 


(12.31 ±0.76 ±0.77) ps"\ 
(8.37 ±0.16 ±0.34) ps^ 1 , 


Cpmp( Lj G eV ) = 1A71 ± Q Q95 ± 0 1Q9 
Ca c pp( 7 GeV 2 ) 


(87) 

( 88 ) 


(89) 


where the first uncertainty is statistical and the second uncertainty is systematic. Together with experimental data, 
Eqs. (87), (88), and (89) will allow determinations of |14fc|, IWftli and |V^,ft/V^ft| with theory uncertainties of 4.4%, 
2.2%, and 4.9%, respectively. A breakdown of the uncertainties into the individual sources, obtained by applying 
Eq. (84) to the various additional form factor fits discussed at the end of Sec. V, is given in Table XII. 


The predicted total decay rates for all possible lepton flavors are 

T(Aft -a p e-u e )/\V ub \ 2 = (25.7 ± 2.6 ± 4.6) ps” 1 (90) 

T(Aft -a p g~^)/\V ub \ 2 = (25.7 ± 2.6 ± 4.6) ps" 1 , (91) 

r(A b -A p T-v»)l \V ub \ 2 = (17.7 ± 1.3 ± 1.6) ps- 1 , (92) 

T(Aft -A A c e~u e )/\V cb \ 2 = (21.5 ± 0.8 ± 1.1) ps" 1 , (93) 

T(Afc -a A c p~^)/\V cb \ 2 = (21.5 ± 0.8 ± 1.1) ps" 1 , (94) 

r(A b -A A c T-v^/\V cb \ 2 = (7.15 ± 0.15 ± 0.27) ps" 1 . (95) 


Motivated by the 7 Z(D^) puzzle [14], we also provide predictions for the following ratios: 


= 0-3318 ± 0.0074 ± 0.0070, (96) 

r(A.£> y A c 6 v 

Vfl \ = 0-3328 ± 0.0074 ± 0.0070. (97) 

T(A 6 -> A c /i-^) v ; 


QED corrections to the decay rates, which may be relevant at this level of precision, have been neglected here. 
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FIG. 14. Predictions for the A& —» pi vt differential decay rates for ^ = e,/r, r in the Standard Model. The inner bands show 
the statistical uncertainty and the outer bands show the total uncertainty, calculated using Eq. (84). 
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FIG. 15. Predictions for the A;, —> A c £ vt, differential decay rates for i = e, fj,, r in the Standard Model. The inner bands show 
the statistical uncertainty and the outer bands show the total uncertainty, calculated using Eq. (84). 
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Cpm* 7 (15 GeV 2 ) 

Ca 0 pc( 7 GeV 2 ) 

C P pp(15GeW) 
Ca„ U i 7(7 GeV 2 ) 

Statistics 

6.2 

1.9 

6.5 

Finite volume 

5.0 

2.5 

4.9 

Continuum extrapolation 

3.0 

1.4 

2.8 

Chiral extrapolation 

2.6 

1.8 

2.6 

RHQ parameters 

1.4 

1.7 

2.3 

Matching & improvement 

1.7 

0.9 

2.1 

Missing isospin breaking/QED 

1.2 

1.4 

2.0 

Scale setting 

1.7 

0.3 

1.8 

z expansion 

1.2 

0.2 

1.3 

Total 

8.8 

4.5 

9.8 


TABLE XII. Approximate breakdown of relative uncertainties (in %) in the partially integrated A& —> ppTu^ and Aj, —¥ 
decay rates and their ratio, defined in Eqs. (87), (88), and (89). As explained in the main text, the combined uncertainty is 
not simply the quadratic sum of the individual uncertainties. 


VII. SUMMARY 

We have presented a high-precision lattice QCD calculation of the complete set of relativistic form factors describing 
the A & —t p and A& —>■ A c matrix elements of the vector and axial vector b -A u and b —» c currents. The form 
factors and their uncertainties in the physical limit are shown in Figs. 10 and 12. Any observable depending on the 
form factors can be calculated using Eq. (84), which is based on two different sets of form factor parameters. The 
“nominal” form factors are used to calculate the central value and statistical uncertainty of the observable, and are 
given by the functions (79) with parameters and correlation matrices from Tables VIII and IX, together with the 
pole masses from Table VII. The “higher order” form factors are additionally needed to calculate the systematic 
uncertainty of the observable, and are given by Eq. (81) with the parameters from Tables X and XI. The higher-order 
fit was performed in such a way that the systematic uncertainty obtained from Eq. (84) includes the continuum 
extrapolation uncertainty, the chiral extrapolation uncertainty, the kinematic ( q 2 ) extrapolation uncertainty, the 
perturbative matching/improvement uncertainty, the uncertainty due to the finite lattice volume, and the uncertainty 
from the missing isospin breaking effects. The individual contributions to the systematic uncertainties in the form 
factors are shown in Figs. 11 and 13. 

Our predictions of the Af, — 1 p t~ and Af, -A A c i~ differential decay rates using the new form factors are 

presented in Sec. VI. The results (87), (88), and (89) for the Af, — > p and At, —► A c differential decay rates 

in the high-g 2 region can be combined with forthcoming experimental data to determine \V u b\, Vcb |, and V u b / U c f, 
with theory uncertainties of 4.4%, 2.2%, and 4.9%, respectively. These uncertainties are competitive with the total 
uncertainties in the 2014 PDG values based on exclusive B meson decays [see Eq. (1)]. Compared to Ref. [29], we 

have reduced the uncertainty in the Af, decay rate at high q 2 by a factor of 3. This reduction in uncertainty 

mainly resulted from the elimination of the static approximation for the b quark. Combined with experimental data, 
our form factor results will also provide novel constraints on right-handed couplings beyond the Standard Model 
[7-10]. The constraints from the baryonic decays nicely complement existing constraints from mesonic decays due to 
the unique dependence of the baryonic decays on ep. Using our Af, -A A c form factors, very precise predictions can 
also be made for the decay Af, -A A c t ~ v t . which may provide new insights into the 7 Z(D^) puzzle [28, 70]. 
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NOTE ADDED 


After the completion of this work, a measurement of the ratio of partially integrated A b —* p/r u tl and A& —> A c /x u tl 
decay rates was published by the LHCb Collaboration, with the result [72] 


f9nax dr(A b ^p H Vp) 

J 15 GeV 2 d q 2 

r?max dTpy-^Ac 

J7GeV 2 d q 2 


d q 2 
d q 2 


(1.00 ±0.04 ±0.08) x 10“ 2 , 


(98) 


where the first uncertainty is statistical and the second uncertainty is systematic. Combined with our lattice QCD 
result in Eq. (89), this gives [72] 


\V ub \ 

\Vc\ 


0.083 ± 0.004(expt) ± 0.004(lattice), 


and, taking the value of \V cb \ extracted from exclusive B decays [72], 

\V ub \ = (3.27 ± 0.15(expt) ± 0.16(lattice) ± 0.06(|y cb |)) x 10" 3 . 


(99) 


( 100 ) 
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Appendix A: Tables of lattice form factor data 


f(A b -»• p) 

|p'| 2 /(2 v/Lf 

C14 

C24 

C54 

F23 

F43 

F63 

f+ 

1 

1.436(60) 

1.417(52) 

1.429(54) 

1.422(60) 

1.419(51) 

1.436(41) 


2 

1.209(66) 

1.202(59) 

1.213(55) 

1.210(83) 

1.218(57) 

1.236(44) 


3 

1.037(73) 

1.037(52) 

1.050(50) 

1.032(56) 

1.051(36) 

1.083(34) 


4 

0.912(64) 

0.925(27) 

0.938(27) 

0.968(31) 

0.964(22) 

0.969(19) 


5 

0.809(34) 

0.823(26) 

0.836(26) 

0.856(30) 

0.857(22) 

0.875(19) 


6 

0.740(34) 

0.754(26) 

0.768(26) 

0.780(29) 

0.788(22) 

0.810(19) 


8 

0.614(35) 

0.651(27) 

0.664(27) 

0.682(29) 

0.690(21) 

0.710(19) 


9 

0.590(36) 

0.608(28) 

0.623(27) 

0.655(29) 

0.649(22) 

0.672(19) 

f± 

1 

1.767(88) 

1.762(64) 

1.802(66) 

1.780(81) 

1.771(73) 

1.804(59) 


2 

1.526(85) 

1.523(68) 

1.558(64) 

1.547(97) 

1.554(71) 

1.582(50) 


3 

1.32(11) 

1.325(85) 

1.353(81) 

1.31(11) 

1.333(76) 

1.390(57) 


4 

1.136(95) 

1.156(61) 

1.185(54) 

1.201(57) 

1.203(34) 

1.224(46) 


5 

1.009(50) 

1.024(39) 

1.056(39) 

1.068(45) 

1.079(32) 

1.112(29) 


6 

0.923(50) 

0.941(39) 

0.968(39) 

0.964(44) 

0.988(32) 

1.029(29) 


8 

0.756(50) 

0.805(40) 

0.830(39) 

0.839(45) 

0.857(32) 

0.893(30) 


9 

0.726(51) 

0.754(41) 

0.781(40) 

0.808(46) 

0.811(33) 

0.849(31) 

fo 

1 

1.056(24) 

1.011(26) 

1.008(26) 

1.051(29) 

1.040(24) 

1.025(21) 


2 

0.887(37) 

0.878(34) 

0.874(37) 

0.889(43) 

0.894(36) 

0.891(36) 


3 

0.775(39) 

0.777(24) 

0.777(24) 

0.788(29) 

0.796(24) 

0.798(23) 


4 

0.718(40) 

0.727(25) 

0.724(29) 

0.757(49) 

0.748(48) 

0.746(20) 


5 

0.654(24) 

0.674(20) 

0.672(20) 

0.666(26) 

0.669(21) 

0.681(17) 


6 

0.606(24) 

0.631(20) 

0.633(19) 

0.626(25) 

0.627(21) 

0.645(17) 


8 

0.516(26) 

0.549(22) 

0.555(20) 

0.550(29) 

0.555(24) 

0.581(19) 


9 

0.502(30) 

0.529(28) 

0.535(20) 

0.530(29) 

0.530(26) 

0.559(21) 

9+ 

1 

0.952(17) 

0.922(17) 

0.920(14) 

0.947(24) 

0.939(23) 

0.925(20) 


2 

0.828(19) 

0.815(17) 

0.813(16) 

0.823(27) 

0.826(24) 

0.821(21) 


3 

0.721(28) 

0.721(19) 

0.722(19) 

0.721(33) 

0.729(25) 

0.735(22) 


4 

0.648(24) 

0.664(16) 

0.663(12) 

0.665(27) 

0.659(24) 

0.663(21) 


5 

0.582(20) 

0.606(15) 

0.606(11) 

0.598(27) 

0.600(27) 

0.607(23) 


6 

0.540(20) 

0.567(15) 

0.569(11) 

0.559(30) 

0.561(31) 

0.576(26) 


8 

0.463(22) 

0.497(16) 

0.505(12) 

0.494(38) 

0.496(36) 

0.517(30) 


9 

0.458(25) 

0.481(24) 

0.489(12) 

0.487(30) 

0.483(30) 

0.507(25) 

9± 

1 

0.952(23) 

0.920(20) 

0.919(16) 

0.947(30) 

0.939(28) 

0.924(25) 


2 

0.827(25) 

0.812(20) 

0.811(18) 

0.818(32) 

0.822(30) 

0.817(26) 


3 

0.719(32) 

0.720(22) 

0.720(22) 

0.715(38) 

0.723(31) 

0.729(27) 


4 

0.643(29) 

0.659(20) 

0.658(16) 

0.657(33) 

0.651(30) 

0.655(26) 


5 

0.578(26) 

0.603(19) 

0.602(14) 

0.586(34) 

0.589(33) 

0.598(28) 


6 

0.535(27) 

0.566(20) 

0.566(15) 

0.546(38) 

0.549(37) 

0.566(32) 


8 

0.456(36) 

0.495(24) 

0.503(17) 

0.476(51) 

0.480(48) 

0.506(40) 


9 

0.454(39) 

0.482(36) 

0.491(17) 

0.473(45) 

0.470(45) 

0.499(36) 

90 

1 

1.475(60) 

1.469(47) 

1.477(41) 

1.496(75) 

1.469(57) 

1.476(43) 


2 

1.237(44) 

1.229(39) 

1.242(36) 

1.274(50) 

1.262(42) 

1.275(26) 


3 

1.055(24) 

1.048(21) 

1.065(21) 

1.082(28) 

1.089(22) 

1.113(17) 


4 

0.912(22) 

0.924(20) 

0.935(19) 

0.982(25) 

0.972(19) 

0.972(15) 


5 

0.808(21) 

0.818(18) 

0.833(19) 

0.872(22) 

0.867(17) 

0.881(14) 


6 

0.734(21) 

0.746(18) 

0.759(19) 

0.782(21) 

0.789(17) 

0.811(14) 


8 

0.603(30) 

0.631(18) 

0.646(19) 

0.668(20) 

0.675(17) 

0.695(14) 


9 

0.576(24) 

0.592(20) 

0.607(20) 

0.640(23) 

0.638(18) 

0.660(16) 


TABLE XIII. At, —> p helicity form factors. 
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/(At -»• p) 

| P '| 2 /(2t r/L) 2 

C14 

C24 

C54 

F23 

F43 

F63 

fY 

1 

1.168(42) 

1.144(46) 

1.152(47) 

1.123(48) 

1.138(37) 

1.164(31) 


2 

0.974(55) 

0.968(54) 

0.977(51) 

0.952(75) 

0.973(49) 

1.000(41) 


3 

0.846(49) 

0.844(31) 

0.858(31) 

0.838(24) 

0.860(20) 

0.889(21) 


4 

0.771(44) 

0.782(19) 

0.792(19) 

0.816(24) 

0.814(25) 

0.819(15) 


5 

0.692(25) 

0.707(20) 

0.715(20) 

0.728(26) 

0.729(21) 

0.745(15) 


6 

0.641(25) 

0.654(19) 

0.666(20) 

0.677(23) 

0.681(20) 

0.698(15) 


8 

0.547(26) 

0.579(20) 

0.590(20) 

0.606(23) 

0.611(19) 

0.629(15) 


9 

0.530(27) 

0.545(21) 

0.557(20) 

0.587(23) 

0.578(20) 

0.597(15) 

fX 

1 

0.505(46) 

0.520(29) 

0.543(25) 

0.556(45) 

0.532(39) 

0.534(30) 


2 

0.465(32) 

0.466(28) 

0.485(25) 

0.503(39) 

0.489(30) 

0.485(23) 


3 

0.397(53) 

0.405(48) 

0.413(43) 

0.397(78) 

0.397(58) 

0.417(32) 


4 

0.307(44) 

0.315(47) 

0.328(46) 

0.325(54) 

0.327(40) 

0.337(40) 


5 

0.267(24) 

0.267(20) 

0.285(18) 

0.288(29) 

0.295(25) 

0.306(18) 


6 

0.238(24) 

0.241(20) 

0.251(18) 

0.243(28) 

0.258(24) 

0.276(18) 


8 

0.176(24) 

0.190(20) 

0.200(18) 

0.198(29) 

0.207(25) 

0.220(19) 


9 

0.165(25) 

0.175(21) 

0.187(19) 

0.187(30) 

0.196(26) 

0.210(20) 

fX 

1 

-0.145(59) 

-0.172(46) 

-0.188(42) 

-0.092(64) 

-0.127(43) 

-0.181(29) 


2 

-0.118(38) 

-0.123(41) 

-0.141(33) 

-0.085(59) 

-0.106(33) 

-0.149(24) 


3 

-0.101(37) 

-0.095(39) 

-0.116(34) 

-0.072(41) 

-0.091(31) 

-0.131(25) 


4 

-0.079(39) 

-0.081(43) 

-0.102(41) 

-0.088(61) 

-0.099(40) 

-0.111(27) 


5 

-0.059(33) 

-0.051(34) 

-0.068(31) 

-0.097(39) 

-0.094(29) 

-0.100(28) 


6 

-0.056(34) 

-0.037(34) 

-0.055(29) 

-0.083(40) 

-0.087(30) 

-0.087(23) 


8 

-0.056(44) 

-0.054(42) 

-0.064(33) 

-0.101(48) 

-0.101(38) 

-0.086(29) 


9 

-0.053(58) 

-0.031(44) 

-0.043(34) 

-0.107(50) 

-0.092(46) 

-0.072(35) 

fX 

1 

0.959(18) 

0.957(19) 

0.950(17) 

0.953(21) 

0.947(17) 

0.953(15) 


2 

0.846(18) 

0.837(20) 

0.838(17) 

0.868(24) 

0.868(22) 

0.865(20) 


3 

0.731(13) 

0.726(12) 

0.738(13) 

0.762(18) 

0.769(16) 

0.773(16) 


4 

0.668(14) 

0.685(15) 

0.687(13) 

0.701(15) 

0.696(13) 

0.701(11) 


5 

0.598(12) 

0.614(12) 

0.621(12) 

0.642(14) 

0.639(12) 

0.643(11) 


6 

0.554(12) 

0.568(11) 

0.576(11) 

0.597(14) 

0.599(15) 

0.606(12) 


8 

0.480(21) 

0.500(11) 

0.508(14) 

0.532(16) 

0.531(13) 

0.543(12) 


9 

0.467(16) 

0.477(13) 

0.487(16) 

0.512(15) 

0.508(14) 

0.522(12) 

fX 

1 

0.008(39) 

0.046(33) 

0.039(29) 

0.007(41) 

0.010(37) 

0.036(34) 


2 

0.024(38) 

0.030(33) 

0.034(34) 

0.061(45) 

0.057(40) 

0.059(41) 


3 

0.015(39) 

0.007(30) 

0.023(33) 

0.058(52) 

0.057(42) 

0.055(41) 


4 

0.030(40) 

0.032(34) 

0.036(29) 

0.053(42) 

0.056(37) 

0.057(37) 


5 

0.024(38) 

0.012(29) 

0.023(27) 

0.068(40) 

0.061(37) 

0.056(34) 


6 

0.023(39) 

0.003(30) 

0.012(27) 

0.062(42) 

0.061(38) 

0.050(35) 


8 

0.030(64) 

0.006(35) 

0.007(31) 

0.068(55) 

0.064(52) 

0.046(42) 


9 

0.017(53) 

-0.006(45) 

-0.004(35) 

0.048(57) 

0.048(57) 

0.028(42) 

fs 

1 

-0.98(12) 

-0.976(89) 

-1.033(80) 

-1.01(15) 

-0.99(12) 

-1.027(86) 


2 

-0.77(10) 

-0.784(98) 

-0.829(94) 

-0.79(12) 

-0.78(11) 

-0.844(79) 


3 

-0.672(44) 

-0.676(39) 

-0.702(32) 

-0.654(53) 

-0.667(43) 

-0.732(33) 


4 

-0.532(43) 

-0.526(37) 

-0.558(31) 

-0.605(53) 

-0.601(42) 

-0.610(32) 


5 

-0.479(40) 

-0.471(35) 

-0.501(29) 

-0.520(48) 

-0.521(40) 

-0.561(31) 


6 

-0.429(40) 

-0.428(34) 

-0.452(29) 

-0.437(45) 

-0.457(40) 

-0.504(31) 


8 

-0.323(41) 

-0.350(36) 

-0.372(30) 

-0.353(47) 

-0.379(42) 

-0.412(33) 


9 

-0.299(43) 

-0.318(37) 

-0.340(31) 

-0.348(50) 

-0.358(44) 

-0.391(35) 


TABLE XIV. A b —> p Weinberg form factors. 
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/(At —> A c ) 

|p'| 2 /(2 tt/L) 2 

C14 

C24 

C54 

F23 

F43 

F63 

/+ 

1 

1.0401(82) 

1.0114(81) 

1.0170(77) 

1.068(15) 

1.062(13) 

1.047(10) 


2 

0.9887(74) 

0.9613(76) 

0.9676(72) 

1.016(14) 

1.011(12) 

0.998(10) 


3 

0.9418(69) 

0.9167(73) 

0.9229(69) 

0.968(13) 

0.963(12) 

0.9528(95) 


4 

0.9018(92) 

0.8826(77) 

0.8852(71) 

0.936(13) 

0.927(12) 

0.9133(88) 


5 

0.864(12) 

0.8515(73) 

0.8528(67) 

0.894(17) 

0.887(17) 

0.876(14) 


6 

0.8327(67) 

0.820(11) 

0.8202(81) 

0.859(15) 

0.852(14) 

0.8455(96) 


8 

0.7692(81) 

0.757(12) 

0.758(10) 

0.799(15) 

0.791(14) 

0.788(11) 


9 

0.7406(98) 

0.727(14) 

0.728(12) 

0.766(15) 

0.758(14) 

0.757(12) 


10 

0.7098(89) 

0.709(16) 

0.708(15) 

0.748(18) 

0.737(16) 

0.735(14) 

f± 

1 

1.467(17) 

1.431(13) 

1.450(17) 

1.464(23) 

1.458(16) 

1.453(14) 


2 

1.400(13) 

1.368(12) 

1.386(13) 

1.398(19) 

1.394(14) 

1.390(13) 


3 

1.339(12) 

1.308(12) 

1.326(13) 

1.335(17) 

1.333(13) 

1.330(12) 


4 

1.268(17) 

1.244(12) 

1.257(15) 

1.282(18) 

1.276(17) 

1.270(17) 


5 

1.219(20) 

1.204(12) 

1.215(12) 

1.228(18) 

1.224(21) 

1.218(21) 


6 

1.180(12) 

1.163(13) 

1.172(12) 

1.185(15) 

1.181(14) 

1.176(11) 


8 

1.094(14) 

1.078(18) 

1.086(15) 

1.103(19) 

1.099(17) 

1.102(16) 


9 

1.056(14) 

1.039(19) 

1.047(16) 

1.055(18) 

1.054(16) 

1.061(16) 


10 

0.997(14) 

1.000(23) 

1.004(21) 

1.021(22) 

1.015(20) 

1.027(20) 

fo 

1 

0.9025(45) 

0.8952(57) 

0.8937(54) 

0.945(13) 

0.9392(98) 

0.9206(62) 


2 

0.8674(41) 

0.8598(54) 

0.8586(51) 

0.906(12) 

0.8996(88) 

0.8846(57) 


3 

0.8336(38) 

0.8273(50) 

0.8258(48) 

0.867(11) 

0.8619(81) 

0.8508(51) 


4 

0.8032(46) 

0.7935(61) 

0.7920(58) 

0.842(13) 

0.8337(86) 

0.8231(75) 


5 

0.7748(41) 

0.7714(46) 

0.7692(43) 

0.805(12) 

0.8012(99) 

0.7939(68) 


6 

0.7458(46) 

0.7429(53) 

0.7409(53) 

0.772(10) 

0.7691(88) 

0.7669(76) 


8 

0.6970(70) 

0.6928(75) 

0.6924(78) 

0.729(17) 

0.723(13) 

0.723(11) 


9 

0.6716(69) 

0.6655(68) 

0.6655(70) 

0.696(14) 

0.691(11) 

0.693(10) 


10 

0.6582(84) 

0.659(12) 

0.658(12) 

0.694(22) 

0.685(18) 

0.685(16) 

9+ 

1 

0.8397(32) 

0.8334(62) 

0.8318(56) 

0.8724(73) 

0.8673(57) 

0.8512(52) 


2 

0.8069(33) 

0.8024(70) 

0.7998(56) 

0.8426(71) 

0.8361(54) 

0.8173(30) 


3 

0.7777(30) 

0.7738(58) 

0.7718(53) 

0.8099(66) 

0.8031(49) 

0.7887(29) 


4 

0.7527(26) 

0.7498(68) 

0.7476(54) 

0.7783(58) 

0.7728(43) 

0.7633(28) 


5 

0.7268(28) 

0.7232(56) 

0.7217(50) 

0.7503(56) 

0.7442(42) 

0.7378(30) 


6 

0.7023(34) 

0.6978(59) 

0.6965(53) 

0.7208(57) 

0.7155(49) 

0.7147(32) 


8 

0.6595(45) 

0.6537(65) 

0.6536(56) 

0.6731(65) 

0.6680(54) 

0.6732(45) 


9 

0.6402(52) 

0.6311(72) 

0.6319(70) 

0.650(11) 

0.6452(91) 

0.6547(62) 


10 

0.6240(71) 

0.624(11) 

0.625(11) 

0.641(10) 

0.6352(87) 

0.642(11) 

9± 

1 

0.8389(35) 

0.8332(50) 

0.8315(48) 

0.8720(72) 

0.8663(57) 

0.8510(54) 


2 

0.8054(37) 

0.8016(67) 

0.7989(52) 

0.8415(72) 

0.8347(57) 

0.8167(34) 


3 

0.7756(35) 

0.7729(51) 

0.7708(44) 

0.8081(68) 

0.8017(50) 

0.7878(33) 


4 

0.7511(32) 

0.7500(64) 

0.7472(55) 

0.7775(56) 

0.7714(42) 

0.7634(32) 


5 

0.7244(34) 

0.7226(47) 

0.7206(48) 

0.7482(54) 

0.7426(41) 

0.7373(33) 


6 

0.6983(41) 

0.6958(50) 

0.6944(45) 

0.7159(56) 

0.7116(47) 

0.7120(33) 


8 

0.6540(57) 

0.6511(57) 

0.6510(49) 

0.6664(56) 

0.6626(48) 

0.6691(41) 


9 

0.6335(58) 

0.6275(67) 

0.6284(56) 

0.6408(94) 

0.6374(81) 

0.6495(52) 


10 

0.6206(90) 

0.6230(91) 

0.6233(95) 

0.637(11) 

0.6325(88) 

0.641(12) 

90 

1 

0.9771(97) 

0.959(12) 

0.9608(99) 

1.007(15) 

0.998(13) 

0.9801(93) 


2 

0.9296(66) 

0.913(11) 

0.9151(88) 

0.958(13) 

0.951(11) 

0.9292(75) 


3 

0.8866(63) 

0.873(11) 

0.8740(86) 

0.916(12) 

0.908(10) 

0.8891(67) 


4 

0.8478(61) 

0.838(11) 

0.8383(86) 

0.873(13) 

0.866(13) 

0.853(11) 


5 

0.8141(62) 

0.804(11) 

0.8044(86) 

0.838(12) 

0.831(10) 

0.8197(88) 


6 

0.7854(66) 

0.774(12) 

0.7737(95) 

0.807(12) 

0.800(11) 

0.7950(76) 


8 

0.7281(86) 

0.716(14) 

0.717(11) 

0.745(15) 

0.738(13) 

0.7409(95) 


9 

0.704(11) 

0.688(16) 

0.690(14) 

0.718(19) 

0.712(16) 

0.717(12) 


10 

0.6732(91) 

0.676(18) 

0.675(16) 

0.698(18) 

0.691(15) 

0.695(13) 


TABLE XV. At —> A c helicity form factors. 
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f(Ab —► A c ) 

Ip'I7( 2 */l) 2 

C14 

C24 

C54 

F23 

F43 

F63 

fY 

i 

0.9534(80) 

0.9273(80) 

0.9316(71) 

0.988(15) 

0.981(14) 

0.965(11) 


2 

0.9096(72) 

0.8842(74) 

0.8894(66) 

0.943(14) 

0.937(13) 

0.923(11) 


3 

0.8697(66) 

0.8465(71) 

0.8517(63) 

0.902(13) 

0.896(12) 

0.8852(95) 


4 

0.8390(82) 

0.8214(75) 

0.8230(65) 

0.877(14) 

0.867(12) 

0.8527(88) 


5 

0.807(11) 

0.7950(70) 

0.7957(60) 

0.840(18) 

0.832(16) 

0.821(13) 


6 

0.7798(63) 

0.768(11) 

0.7679(80) 

0.810(16) 

0.802(14) 

0.7954(98) 


8 

0.7255(78) 

0.714(12) 

0.7154(98) 

0.758(15) 

0.749(15) 

0.746(11) 


9 

0.7009(96) 

0.688(13) 

0.689(11) 

0.730(15) 

0.721(14) 

0.719(12) 


10 

0.6761(87) 

0.676(15) 

0.674(14) 

0.716(19) 

0.704(16) 

0.700(14) 

f¥ 

1 

0.365(11) 

0.3569(91) 

0.367(11) 

0.338(18) 

0.338(13) 

0.3455(96) 


2 

0.3481(80) 

0.3427(62) 

0.3516(72) 

0.323(13) 

0.3241(87) 

0.3308(73) 


3 

0.3330(67) 

0.3273(56) 

0.3354(60) 

0.308(11) 

0.3107(83) 

0.3149(71) 


4 

0.3045(82) 

0.2997(70) 

0.3072(83) 

0.287(11) 

0.2903(86) 

0.296(11) 


5 

0.2927(78) 

0.2903(52) 

0.2968(56) 

0.2756(10) 

0.2784(76) 

0.2813(75) 


6 

0.2839(60) 

0.2802(51) 

0.2861(55) 

0.2665(98) 

0.2690(74) 

0.2700(63) 


8 

0.2618(76) 

0.2578(66) 

0.2624(65) 

0.244(11) 

0.248(11) 

0.2523(94) 


9 

0.2524(77) 

0.2489(66) 

0.2531(66) 

0.231(11) 

0.2363(10) 

0.2422(90) 


10 

0.2279(75) 

0.2300(88) 

0.2330(78) 

0.216(13) 

0.221(11) 

0.231(12) 

fY 

1 

-0.090(15) 

-0.057(13) 

-0.068(11) 

-0.075(32) 

-0.074(29) 

-0.079(20) 


2 

-0.078(14) 

-0.045(12) 

-0.057(10) 

-0.070(31) 

-0.070(28) 

-0.072(20) 


3 

-0.070(13) 

-0.037(12) 

-0.051(10) 

-0.067(31) 

-0.065(29) 

-0.067(20) 


4 

-0.073(14) 

-0.057(13) 

-0.064(11) 

-0.072(31) 

-0.067(28) 

-0.060(19) 


5 

-0.068(22) 

-0.051(12) 

-0.0572(97) 

-0.074(41) 

-0.065(40) 

-0.058(31) 


6 

-0.076(14) 

-0.057(17) 

-0.061(11) 

-0.084(31) 

-0.073(28) 

-0.064(21) 


8 

-0.072(17) 

-0.054(21) 

-0.058(14) 

-0.075(50) 

-0.066(44) 

-0.057(28) 


9 

-0.078(25) 

-0.060(28) 

-0.063(21) 

-0.091(44) 

-0.080(38) 

-0.067(27) 


10 

-0.051(17) 

-0.049(20) 

-0.048(15) 

-0.063(54) 

-0.053(46) 

-0.042(29) 

A A 

1 

0.8581(72) 

0.838(12) 

0.8386(87) 

0.881(12) 

0.890(10) 

0.8552(65) 


2 

0.8235(57) 

0.810(12) 

0.8089(85) 

0.853(10) 

0.8514(86) 

0.8242(58) 


3 

0.7918(53) 

0.780(11) 

0.7791(84) 

0.822(10) 

0.8126(84) 

0.7948(58) 


4 

0.7608(51) 

0.749(11) 

0.7498(83) 

0.7824(94) 

0.7795(83) 

0.7627(69) 


5 

0.7357(55) 

0.726(11) 

0.7255(84) 

0.7579(96) 

0.7502(81) 

0.7397(67) 


6 

0.7144(56) 

0.704(12) 

0.7030(96) 

0.7355(95) 

0.7271(82) 

0.7232(68) 


8 

0.6709(71) 

0.659(13) 

0.659(10) 

0.687(12) 

0.6791(98) 

0.6818(88) 


9 

0.6519(90) 

0.637(15) 

0.638(14) 

0.666(17) 

0.659(14) 

0.664(11) 


10 

0.6289(80) 

0.626(16) 

0.627(14) 

0.647(14) 

0.639(11) 

0.643(12) 

fY 

1 

0.032(12) 

0.008(15) 

0.012(10) 

0.016(18) 

0.040(15) 

0.007(14) 


2 

0.031(11) 

0.014(15) 

0.017(10) 

0.020(15) 

0.028(13) 

0.013(12) 


3 

0.027(11) 

0.013(15) 

0.014(10) 

0.024(15) 

0.018(12) 

0.012(11) 


4 

0.016(10) 

-0.002(14) 

0.0045(98) 

0.008(16) 

0.014(14) 

-0.001(15) 


5 

0.019(11) 

0.005(15) 

0.008(10) 

0.016(15) 

0.013(13) 

0.004(14) 


6 

0.027(12) 

0.014(18) 

0.015(13) 

0.033(15) 

0.026(13) 

0.019(13) 


8 

0.029(15) 

0.013(24) 

0.014(17) 

0.034(18) 

0.028(16) 

0.022(15) 


9 

0.031(17) 

0.017(27) 

0.016(21) 

0.043(22) 

0.037(20) 

0.025(18) 


10 

0.014(16) 

0.006(20) 

0.005(13) 

0.016(20) 

0.011(17) 

0.003(15) 

fY 

1 

-0.501(24) 

-0.513(46) 

-0.525(43) 

-0.530(42) 

-0.457(34) 

-0.531(30) 


2 

-0.467(16) 

-0.459(14) 

-0.478(14) 

-0.462(29) 

-0.438(24) 

-0.467(16) 


3 

-0.438(15) 

-0.431(13) 

-0.447(12) 

-0.433(27) 

-0.439(23) 

-0.439(15) 


4 

-0.422(17) 

-0.437(17) 

-0.438(15) 

-0.439(37) 

-0.421(32) 

-0.440(26) 


5 

-0.399(16) 

-0.402(14) 

-0.410(13) 

-0.408(28) 

-0.411(25) 

-0.410(18) 


6 

-0.381(16) 

-0.378(14) 

-0.387(13) 

-0.384(27) 

-0.388(24) 

-0.387(16) 


8 

-0.343(26) 

-0.342(15) 

-0.352(14) 

-0.348(33) 

-0.348(30) 

-0.354(17) 


9 

-0.330(31) 

-0.329(17) 

-0.337(16) 

-0.331(29) 

-0.331(26) 

-0.339(19) 


10 

-0.299(31) 

-0.339(35) 

-0.332(31) 

-0.348(59) 

-0.345(52) 

-0.355(43) 


TABLE XVI. Ab —> A c Weinberg form factors. 
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